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Abstract 

We discuss the choice of weight functions for the moments of the spectral density for two- 
point correlators of hadronic currents over a finite energy interval. Of phenomenological 
relevance is an analysis of the spectra of r lepton decay on the energy interval [0, M T ] and 
the low-energy hadron production in e + e~ annihilation. General arguments are given for 
the calculability of such moments in perturbation theory using both a finite-order analysis 
and infinite resummation on the contour. Nonperturbative contributions emerging from 
the operator product expansion for two-point correlators are discussed within explicit 
models for the physical spectra. The quantitative analysis strongly disfavours weight 
functions that suppress the high-energy contribution to theoretical moments. This is in 
agreement with expectations from qualitative considerations in perturbative QCD with 
asymptotic freedom. We discuss the implication of our results for the ultimate accuracy 
that can be reached in r decays and low-energy e + e~ annihilation into hadrons with 
present experimental data. 
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1 Introduction 



The central quantities of interest for studying strong interactions in low-energy hadron phe- 
nomenology at accelerators are the spectra of produced particles as being related to the cross 
section and described by spectral functions. The modern theory of strong interactions (QCD) 
describes the production processes within perturbative expansions in the strong coupling con- 
stant in terms of quarks and gluons and is unable to predict the low-energy hadron spectra 
pointwise as functions of the energy. Only some integral characteristics of low-energy hadron 
spectra such as the average over an energy interval can be reliably calculated within pertur- 
bative QCD. Theoretically, the procedure of averaging is best justified for observables related 
to two-point correlators of hadronic currents with well established and simple analytic proper- 
ties as functions of the energy. Phenomenologically, this is the case for e + e~ annihilation into 
hadrons and hadronic r lepton decays where experimental data of high precision are available. 

The process of e + e~ annihilation into hadrons has been a laboratory for studying strong 
interactions since long ago where important results have been obtained (as a recent review, see 
e.g. Ref. As a famous example of such results one can mention the constraint on the number 
of colours N c in QCD obtained from the normalized decay rate R e + e - (s) of e + e~ annihilation 
into hadrons smoothed over a finite interval at low energies of about y/s ~ 2 4- 3 GeV 



where a s is a strong coupling constant. During the last decade much attention has also been 
paid to a study of hadronic (semileptonic) r lepton decays which provides a wealth of informa- 
tion on low-energy hadronic physics where the accuracy of experimental data is permanently 
improving || |3|. This makes r lepton physics an important area of high-precision particle 
phenomenology. 

In the present paper we discuss in some detail the choice of the moments which are given 
by integrals of the spectral density over the energy range available in experiments. These 
integrals may contain specific weight functions of different shape that determine the particular 
form of averaging and the kind of moments. The moments with suitable weight functions 
can be computed theoretically in perturbative QCD while the detailed pointwise description 
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of low-energy spectral functions themselves as functions of the energy requires a quantitative 
understanding of the hadronization mechanism which is lacking at present, especially for light 
hadrons. The procedure of averaging, in its different modifications, is referred to as a property 
of duality between hadronic and quark-gluon spectra j4|. 

While generating a particular set of the moments of the spectral density as observables 
to confronting theory with experiment, one faces the usual conflict between the precision re- 
quirements dictated by theory and experiment. In perturbative QCD a reliable computation 
is possible only at large energies where the perturbation series converges in a sense that a fi- 
nite number of high-order corrections in the strong coupling constant (usually one or two) is 
considered to be numerically small because of the property of asymptotic freedom, so that the 
moments which emphasize the high-energy part of the spectrum are most reliably calculated. 
However, such moments are not welcome from the experimental point of view because contri- 
butions coming from the high-energy part of the spectrum in the range < s < usually 
have got a poor experimental accuracy. On the other hand, moments which suppress the high- 
energy part of the spectrum are saturated by low-lying resonances or few-particle states which 
are easier to measure experimentally compared to the states with many particles. This gives 
rise to better experimental precision for such moments, which would induce on using these 
moments for a theoretical analysis. However, the theoretical uncertainty for such moments is 
larger because of both the poor convergence of the perturbation series due to the large value of 
the effective coupling and the large nonperturbative contributions stemming from the infrared 
region which are difficult to estimate quantitatively. 

The semiphenomenological way for estimating the nonperturbative contributions based on 
the operator product expansion (OPE) for two-point correlators is most advanced where the 
knowledge of the numerical values for a few low-dimension terms (gluon condensate and four- 
quark condensates) fixed from the data. Higher order terms of dimensionality larger than six are 
virtually unknown; the numerical values that are commonly used derive from some approxima- 
tions (mainly the factorization into condensates of lower dimensionality) or models (instantons 
for gluonic condensates) with theoretical uncertainties being scattered over a wide range with 
rather hypothetical justification in the literature. This makes the estimate of contributions of 
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nonperturbative terms of higher orders rather imprecise. 

Nowadays it is a general trend to consider restrictions imposed by perturbation theory con- 
vergence to become more important than those coming from nonperturbative terms. With the 
knowledge of higher and higher order terms of the perturbation series the asymptotic character 
of perturbation theory is more and more pronounced which forces one to move the renormal- 
ization point to larger and larger energies. In this region nonperturbative terms die out fast 
because being power corrections they decrease faster than the terms of the logarithmic pertur- 
bation series. On the other hand, power corrections reveal the general structure of observables 
computed in QCD, i.e. they account for nonperturbative effects. In this circumstances there is a 
tendency to introduce power corrections almost by brute force in case where OPE is not directly 
applicable because their existence seems to be justified by general principles. Even though the 
motivation for such a development is clear, the present state of the analysis is definitely far 
from being quantitative. The obstacle is also that this kind of analysis is infrared dependent, 
i.e. it depends on the properties of strong interactions at small energies where perturbation 
theory - the only universal tool of quantitative investigation - is not applicable. This fact 
makes conclusions based on power corrections depend on the particular model or convention 
(for instance on a recipe of resummation). 

The main purpose of this paper is to give explicit arguments and quantitative estimates for 
the calculability of the moments of two-point correlators in perturbative QCD and to find a 
compromise choice of weight functions for the moments that optimizes the precision comparison 
of theoretical input with experimental data. To be specific, we are talking about the precision 
analysis of r lepton decays having been of much interest recently where experimental data are 
very good. Therefore, we concentrate on the aspects of accuracy of theoretical calculations for 
the QCD part of the differential r lepton decay rate and its moments. It is the high precision 
achieved in the experimental analysis of r decays and the rather advanced stage of theoretical 
description that calls for a critical examination of the theoretical formulae used in the analysis 
of the physics of the r system. In particular the criteria of reliability of theoretical formulae 
should be considered carefully. The quantitative discussion of this issue is possible because 
of the simplicity of the Green's functions necessary to describe the process - only two-point 
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correlators depending on one single energy variable are relevant. The analytic properties of the 
two-point correlators in the energy variable are strictly fixed from general principles of quantum 
field theory. An additional simplification for the renormalization group analysis is provided by 
the fact that only one scale is involved in the theoretical description. 

Theoretical calculations for the r system can be done within OPE for two-point correlators 
which contain perturbative and power corrections. Perturbation theory can then be further 
improved using methods of summation of an infinite number of terms within the renormalization 
group. We consider these steps in turn. In Sec. 2 we give the basic ideas and fix the notation. In 
Sec. 3 we consider finite-order perturbation theory, in Sec. 4 we consider power corrections. As 
an important phenomenological example we consider some peculiarities related to quark mass 
corrections to current correlators which are important for the determination of the strange 
quark mass (Sec. 5). In Sec. 6 we consider the infinite resummation of perturbative terms 
generated by the running of the coupling constant within a renormalization group approach. 
Sec. 7 contains our conclusions. 



2 Setting the stage: details of the theoretical description 

Having described the motivation of our work which is rather general and valid for two-point 
correlators of hadronic currents, we now give the specific details for the description of r decay 
observables. Similarities to the e + e~ annihilation into hadrons are obvious and applications to 
other hadronic channels are straightforward. In this section we mainly fix our notation; the 
consideration is rather standard and can be found in the numerous literature on this subject [[| 

I- 

The semileptonic (hadronic) r lepton decay is mediated by the charged weak hadronic 
current of the form 

j"(x) = Vudtij^l - 7 5 )rf + VusU'j^l - 7 5 )s (1) 

where V u d and V us are Cabibbo-Kobayashi-Maskawa matrix elements (elements of the weak 
mixing matrix). The correlator for the weak hadronic currents in Eq. ([I]) has the general form 

IV(g 2 ) = 12tt 2 z J (T JlM (x)ji(0))e^dx = q^U q (q 2 ) + g, u Il g (q 2 ) (2) 
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where U q (q 2 ) and U g (q 2 ) are invariant scalar functions. These scalar functions are further 
specified depending on which current is considered. In case of the (ud) quark current j^x) = 
«7 M (1 — 75)6/ (denoted as light quark case) the massless limit is assumed in which case the 
correlator is transverse, i.e. both invariant functions U qjg (q 2 ) are expressible through a single 
scalar correlator function U u d (q 2 ) 

U q (q 2 ) = Il ud (q 2 ), U g (q 2 ) = -q 2 Yi ud (q 2 ). (3) 

The correlator in case of the (us) quark current (the term proportional to V us , also referred 
to as the strange quark case) is slightly different as the nonvanishing strange quark mass is 
taken into account. The relevant formulae are given e.g. in Refs. || |7|, [|. We will discuss the 
strange case later. The rest of the consideration is rather similar to the light quark case. We 
skip the specification (ud) in the following and use U u d(q 2 ) = n(g 2 ). We mention that we work 
within QCD with three light quarks and do not consider corrections due to heavy quarks (c 
quark) which would enter the calculation at high orders of perturbation theory through loop 
effects [0], p]] . In addition to the different current specifications as in Eqs. ([I]) and @ it is 
convenient to consider the vector and axial parts of the correlator separately, 

U v+A (q 2 ) = U v (q 2 ) + U A (q 2 ) (4) 

with U v (q 2 ) beging related to the pure vector part and II^(g 2 ) related to the pure axial part. 
We introduce the spectral density of a correlator as the discontinuity across the physical cut, 

1 1 

p(s) = — DiscII(s) = — (U(s + iO) -Il(s-iO)), s>0 (5) 
Z7T2 2m 

which splits into vector and axial-vector parts accordingly, 

Pv+a(s) = pv(s) + Pa(s). (6) 

Here we mainly concentrate on the massless limit for the correlator. The case of the strange 
quark with nonvanishing mass will be discussed later. In the following the indices V and A 
will be omitted in the generic case as well if no confusion arises. The correlator in Eq. (Q) is 
normalized to the number of colours N c (the same holds true for the spectral density which 
means that p(s) — > N c for s — > 00) in the leading parton model approximation with massless 
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quarks. In the following we will occasionally use also a slightly different normalization which 
explicitly accounts for the number of colours, resulting in a correlator which is normalized to 
unity. 

The above considerations are general and are also used by experimentalists to classify the 
appropiate channels: strange particles (K mesons) form the strange channel, non-strange axial- 
vector mesons (as a\ and the like) form the axial- vector channel, and the classical vector meson 
p represents the non-strange vector channel. In this respect the pion is somewhat special. It is 
a Goldstone boson with spin zero and gives a contribution to the axial correlator IT^ (q 2 ) in the 
massless limit. 

Before specifying the theoretical calculations we give the general form of an important 
observable of r decays. The basic observable is the normalized r lepton decay rate for the 
decay of the r lepton into hadrons written in the standard form as 

1 (r — > t + Vi + V T ) v ' 

The leading terms in Eq. (0) are the parton model results while the terms S u d and S us represent 
the effects of QCD interactions and mass effects (in case of nonvanishing quark masses) |12|, [13 



14], 15J, [L6| . V u d and V us are matrix elements of the weak mixing matrix as defined in Eq. ([[]) , 



and Sew describes the electroweak radiative corrections to the r decay rate fll7| . Now we take 
a deeper look into the physics of hadronic r decays. In the massless limit the expression for 
the t decay rate relevant for the QCD part of the quark current is given by the phase-space 
integral 

B? » = WeJ J- ! (l-^) , ( 1 + 2^)^ W *. (8, 
The spectral density p(s) is related to Adler's function D(Q 2 ) through the dispersion relation 

where Q 2 = —q 2 . At this point theory enters. It is Adler's function D(Q 2 ) that is most 
convenient to compute theoretically in the Euclidean domain. It allows one to theoretically 
predict r observables through a theoretically computable spectral density p{s). Still, to extract 
the theoretical prediction for p(s) from D(Q 2 ) is not straightforward. The point is that D(Q 2 ) 
is only known as an asymptotic expansion at large Euclidean values Q 2 while p(s) is obtained 
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Figure 1: experimentally measured values for the relative cross section in the e + e channel, 



the data set is taken from Ref. [181. 



as a discontinuity of the correlator across the physical cut. Therefore, an analytic continuation 
into the complex g 2 -plane is necessary. Analytic continuation is an improperly posed problem, 
i.e. small errors in the initial function D(Q 2 ) can produce big errors in p(s). This instability 
is especially important for a theoretical computation of p(s) at low energies. This situation 
can also be reformulated in the language of integral equations. The dispersion relation in 
Eq. (H) gives the spectral density p(s) through Adler's function D(Q 2 ) as a solution of the 
integral equation. The integral equation given by Eq. @ is a Fredholm integral equation of 
the second kind which is known to lead to an improperly posed problem. Thus, errors for the 
function p(s) (as solution of this integral equation) are not continuously related to errors for 
D(Q 2 ) (as initial data of the integral equation) and can be very large. The general procedure 
of constructing approximate solutions for such a problem was suggested by Tikhonov and is 
known as smoothing. Averaging the spectral density over a finite energy interval (as in sum 
rules that correspond to the duality concept) can be considered as a particular realization of 
Tikhonov's smoothing procedure. 

The function p(s) at low energies is the main object entering Eq. (§), therefore we want to 
concentrate our studies on it. The counterpart of p(s) on the experimental side is the hadronic 
spectral density p had (s) which can be measured in r decays in the finite squared energy interval 
[0,M?] with M T = 1.777 GeV §. The hadronic spectral density p had (s) is a rapidly varying 
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function in the vicinity of resonances (for the experimental data set in the vector channel see 
e.g. Fig. |I]), therefore there is no hope that one is able to compute this function theoretically 
pointwise within perturbation theory. This fact is related to the so-called confinement (or 
hadronization) problem. Instead of a pointwise description of the spectrum at low energies, the 
appropriate quantities to be analyzed theoretically in perturbative QCD are the moments or 
integrals of the spectrum with a complete set of weight functions. 
We define moments of the spectral density by 

rM% J rig 

m = (j + i)/ o p(*)jM*F* =:1+mi ' (10) 

The factor (Z+l) in the definition of the moments is chosen to have all contribution of the parton 
part uniformly normalized to unity, as it is written explicitely in Eq. (|T0|). This corresponds to 
the parton contribution which is independent of QCD interactions. Equivalently one can say 
that all measures 

defined on the interval [0, M%] are normalized to 1 for the volume of the integration space 
which in this case is the interval [0, M 2 ] ■ The quantities mi defined in Eq. ( |T0D are related to 
the interaction. They are of real interest because they contain a wealth of information about 
the hadronic spectral density p had (s) and are sufficient for theoretical studies. 

Theoretical calculations of the moments can be done within OPE which contain perturbative 
and power corrections. The perturbation theory can then be further improved using methods 
of summation. We consider these steps in turn. 

3 Finite-order perturbation theory 

In this section we start the analysis of the moments given in Eq. ( |I0|) with the consideration 
of the ordinary perturbative part of the theoretical spectrum or Adler's function D(Q 2 ). The 
theoretical prediction for the function D(Q 2 ) has been calculated with a very high degree of 
accuracy within perturbation theory (see e.g. Refs. |TP|, [2D], pifl). 

In the massless limit for light quarks (u, d) the perturbation theory expression for Adler's 
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function is the same for both the axial and the vector channel. In the MS-scheme it reads 



Dm = 1+ ^ +kl ^y +k2 ^y +h ^y + o(aa (i2) 

where the running coupling a s is normalized at the scale /i 2 = Q 2 , a s = a s (Q 2 ). Here 

299 

h = — -9C(3)« 1.63982, 

= 58057 _ 779 75 

288 4 w 2 w v ; 

for = 3 light flavours, ((z) is Riemann's ( function. The fourth-order MS-scheme coefficient 
&3 in Eq. fll2"|) is related to the divergent parts of five-loop diagrams and is not yet known 
at present. Eqs. (|T2|) and (0) constitute the set of theoretical information necessary for a 
perturbation theory analysis of the r system. The analysis of finite-order perturbation theory 



was presented in Ref. |22 



Starting with Adler's function given in Eq. (|T^) we deduce the spectral density p(s) within 
the MS-scheme, 

M = i + 2£1 + h + (, 2 - f^)' + . . . (u) 

where the term proportional to it 2 is a result of the analytic continuation from the Euclidean 
domain, /3 is the leading coefficient of the /3-function. 

The first step in studying the perturbative part of the moments mi in Eq. (|1(]) is to fix our 
criteria on how to decide which weight functions are theoretically preferable. In ordinary (or fi- 
nite order) perturbation theory a natural criterion is the explicit convergence of the perturbation 
series, i.e. the pattern of numerical behaviour of consecutive terms of the perturbation series. In 
general, one expects the convergence of the perturbation series to be better for moments which 
are more 'perturbative', i.e. sensitive to larger scales. However, the explicit convergence can 
be concealed by the use of a particular scheme, for instance by defining the coupling constant 
in the MS-scheme which is most widely used. Therefore, the criterion of explicit convergence 
is scheme dependent. Besides physical reasons of the applicability of perturbation theory to a 
given observable the explicit convergence for the moments (numerical structure of the pertur- 
bation theory series with finite number of terms) can look better or worse depending on the 
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choice of the expansion parameter in a given scheme. This is an artefact of the definition of 
the coupling (renormalization or subtraction). 

To eliminate such an artificial influence of the scheme definition on physical conclusions one 
should use scheme independent criteria of theoretical calculability within perturbation theory. 
The invariant content of the investigation of the spectrum in perturbative QCD, independent 
of any definition of the charge or the coupling parameter, is the simultaneous analysis of a set 
of moments. Technically, an efficient way to do this is to factor out all unrelated constants by 
redefining the charge and introducing an effective coupling which absorbs all constant terms of 
the perturbative expansion for the polarization function or the spectral density. The purpose of 
the effective coupling is to get rid of artificial scheme- dependent constants in the perturbation 
theory expressions. Therefore, in order to analyze a set of moments mi in the most transparent 



way (see e.g. Refs. |23|, g4[ p5| ) we express the spectral density in terms of an effective coupling 
for the t system. 

Note that the use of an effective coupling is a technical trick which simplifies the analysis. 
Such a parameter is necessary in order to be able to assess the explicit convergence of the 
perturbation series. The really invariant measure is given by the mutual relations between 
moments while the effective coupling is still a parameter of the theory to which one cannot 
prescribe any special physical meaning. Despite this fact, the introduction of a natural inter- 
nal coupling parameter allows one to extend the perturbation theory series available for the 
description of the relations between observables by one more term as compared to the analysis 



in e.g. the MS-scheme (see e.g. Refs. [^3|, 



In this section we work within finite-order perturbation theory (FOPT) and define an effec- 
tive coupling ojif(s) on the physical cut for sufficiently large values of s by the relation 

aM(a) = «£> + h (-My + ( fc _ (^) 3 + . . . (is, 

such that 

p(s) = l + a M {s). (16) 

The subscript "M" stands for a Minkowskian definition of the effective coupling, i.e. the defini- 
tion on the physical cut. The decomposition of the spectral density in Eq. ([UJ) reflects the fact 
that within perturbation theory the correlator contains the parton part which is independent 
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of a s . All the constants that may appear in the perturbation theory expression for the spectral 
density p(s) due to a particular choice of the renormalization scheme are absorbed into the 



definition of the effective charge (see e.g. Refs. |2^, |29], |30|), so that only effects of the 
running of the coupling itself are left. The solution of the evolution equation for the effective 
coupling, 

da M (s) 



ds 



P(a M (s)) = -a M (s) (f3 a M {s) + fta M (s) 2 + (3 2 a M {s) 3 + 0(a A/ (s) 4 )) (17) 



resulting from the renormalization group analysis of the correlator with a given effective j3- 
function /3(a) can be obtained by quadrature, 

f^ds' 



an 



s) = a M + [ /3(a M (s'))-^- 
Jui s' 



= a M + /3 La A/ + (/3 1 L + $L> Af + (^ (18) 

where aM = clm{M 2 ) and L = ln(M^/s). Moments of the spectral density with different weight 
functions emphasize different regions of the evolution trajectory of the renormalization group 
equation. Because we are interested in the running of the effective coupling, we can also use this 
solution in terms of an expansion in L which is more relevant for our particular setup, instead 
of an expansion in aM (which is standard). At any given order of perturbation theory for the 
effective or for the exact function (3{a) the running coupling ciAf(s) is given by the evolution 
translation of the initial value aM according to the renormalization group equation. The result 
reads 



a M {s) = a M - p{a M )L + -fj{a M )- 



' da 



6 da da 

a=a M 



L 3 + 0(L 4 ) 

a=a,M 



exp (-Lfta)£)a 



(19) 

a=a,M 

where the solution of Eq. ( TiT| ) is written in a symbolic operator form. Note that this expression 
is useful for the so-called 't Hooft scheme where, by definition, the /3-function of perturbation 
theory is given by the concise expression (3{a) = — j3o(a/n) 2 — /5i(a/7r) 3 pl| . 

Defining the effective coupling gim(s) directly through the spectrum p(s) itself one obtains 
perturbative corrections to the moments in Eq. ( |I0D only because of running. Without running 
(as, for instance, in the conformal limit of QCD with a vanishing /3-function or at the infrared 
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fixed point) one would have 

M z = l + a M (M T 2 ) or = a M {M 2 T ) (20) 

with mi defined in Eq. ([I0|). In this situation the whole theoretical description of the physics of 
the r system in the massless approximation (i.e. without strange particles and including only 
perturbative corrections without possible power corrections) would reduce to the determination 
of a single number a^(M^) and consequently there would be no problems with the convergence 
of the perturbation theory series. This reminds one of the situation in QED where the coupling 
is defined through the subtraction on the mass shell and simply gives the cross section in the 
Thomson limit. But because of the running of om(s) (which is important numerically because 
both ajw(M^) and (3q are large), different observables, i.e. different moments of the spectral 
density, generate different perturbation series from the original object p(s) in Eq. fll6|) and 
this difference is now within the reach of experiments. Thus, moments just allow one to study 
the evolution or /3-function of the effective coupling cim{s). We remind the reader that in 
this section we discuss only the perturbative part of the theoretical spectrum where power 
corrections are neglected. 

The contributions of powers of logarithms (from Eqs. ([18]) or ((L^)) to the normalized mo- 
ments in Eq. (|II]) are given by 

and 

n ^ f M K 2( M r\ S l ds 2 

A general formula for an arbitrary (integer) power of the logarithm reads 

(/+1 U ln (vJ(M^ = (7Tiy- (23) 

Therefore, at any fixed order of the perturbation series expansion the effects of running die out 
for large values of I improving the (explicit) asymptotic structure of the perturbative series for 
the moments in Eq. (p!U|). This is obvious and expected in QCD with its property of asymptotic 
freedom because the moments with weight functions s l suppress the infrared (small s) region 
of integration where perturbation theory is not applicable. Note that for a given order of 
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perturbation theory the real effects of running are connected with powers of logarithms and 
not with powers of the coupling a u- Therefore one computes an expansion of the spectral 
density p(s) as a Taylor series with the initial value at M T , 

p(s) = po(M T ) + Pl {M T )L + p 2 (M T )L 2 + ... (24) 

It is just the coefficients of this Taylor series (coefficients pi(M T ) of the powers of logarithms) 
that are calculated by a perturbative expansion in the coupling constant. In the effective scheme 
using aw one finds 

po(M T ) = a M , 
Pl {M T ) = -P(a M ), 

P2(M T ) = ^(a M )^, ... (25) 
while for the MS-scheme one would have 

MMr) = ^MA + kl {^MA)\ ( k2 - ^ij (^J + . . . (26, 

The technical advantage of introducing the effective scheme becomes obvious. The effects of 
running as given by powers of logarithms can be obtained in a concise form as explicit functions 
of aM from Eq. (|23| ) if the /^-function is known. Note that only for the leading order /3-function 
the powers of logarithms coincide with the powers of the coupling a. 

The behaviour at fixed I and large p, i.e. the convergence (or asymptotic structure) of the 
perturbation series with an infinite number of terms for the moments is another story. This 
question was intensively discussed in the literature. We present our analysis of this problem in 
Sec. H 

Therefore, the set of moments in Eq. ( |lOl) has a simple and clear structure of perturbative 
convergence in finite-order perturbation theory: larger values of / are better behaved from 
the perturbation theory point of view. However, in practice and for an efficient comparison 
with experimental data the choice of weight functions is also dictated by the precision of the 
available experimental data. With this constraint taken into account, s'-moments for large 
values of / are not welcome from an experimental point of view. They are dominated by the 
contributions coming from the high-energy part of the r decay spectrum (and therefore show 
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a better perturbative convergence) while the experimental accuracy for the moments basically 
deteriorates with increasing / because poorly known contributions close to the high-energy end 
of the interval are enhanced. Therefore, with respect to the precision aimed for the the structure 
of the set of the s'-moments on the experimental side is opposite to that on the theory side: 
for larger / the precision for the calculation of moments using experimental data is worsening 
fast. 

To balance the precision requirements from the experimental and theoretical sides one can 
consider a modification of the set of moments (see e.g. Refs. |]32], |33fl). Note that in general any 
polynomial in s can be used as a weight function for a finite interval. In a sense a complete 
system of orthogonal polynomials (like the Jakobian polynomials p( a,h '(x) with the rather 
general weight function x a (l — x) b which are orthogonal on the finite interval (0, 1) ||34|| ) could 
provide an exact expansion of the function p(s) in term of its moments and allows for its full 
reconstruction of the special shape of the coefficients of the expansion (supposing that p(s) is 
a continuous function). Therefore, the particular choice of weight functions is determined by 
the efficiency of solving concrete problems. Actually, the basis {s l ; I = 0, . . . , 00} is complete 
for continuous functions on the interval [0,Mj?], and the choice of this basis allows for the 
most direct evaluation of the structure of the moments from the perturbation theory point of 
view which is the most relevant for the theoretical analysis. At s = the spectrum is quite 
perturbative and for large values of I the s z -moments are basically determined by p(M^). 

As the simplest modification done in order to suppress experimental errors from the high- 
energy end of the spectrum, we introduce the system of modified moments with weight functions 



(fc + Z + l)! / _s_ 

kill [ Ml) \M1 



™«(«) = - '/„; — U-T79 ) ■ (27) 



Even in the case / = we can consider the modified moments as mixed moments, because they 
are just linear combinations of the direct moments with weight functions s l as given in Eq. (III). 
The modified moments 

(fc + l + 1)! f M ? ( __ _s_\ k ( s V p(s)ds 



are normalized to unity. Within the set given in Eq. (p8|) the best choice from the experimental 
point of view is to use large values for k and small values for /. This choice was also advocated 
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to be justified theoretically as improving the precision based on the integration over the contour 
in the complex q ,2 -plane (see Fig. The reasoning was that the weight functions (1 — s/M 2 ) k 
suppress the contribution of that part of the contour that is close to the real positive semi-axis 
where OPE is not applicable (region A in Fig. |^) which in turn, according to standard wisdom, 
can improve the accuracy of theoretical predictions. Below we give arguments that this is 
not the case: the weight functions (1 — s/M 2 ) k simply ruin the perturbative structure of the 
moments for large values of k. Let us demonstrate this in more detail. 

Within finite-order perturbation theory the integration over the contour |35|, |36|, |37 



is com- 



pletely equivalent to the integration over the cut along the positive semi-axis because of the 
analytic properties of the functions ln p (— M 2 /q 2 ) occuring in the expansion of the correlator. 
This can easily be seen from the spectral density according to Eq. (0). All polynomial mo- 
ments can be rewritten as contour integrals in the complex g 2 -plane. Therefore, for finite-order 
perturbation theory it makes no difference which particular representation for the moments is 
used: integrals over the contour in the complex g 2 -plane and integrals over the cut along the 
positive semi-axis are mathematically equivalent at any order of perturbation theory. Power 
corrections and infinite summation in perturbation theory will be discussed later on. 

In the framework of finite-order perturbation theory we will discuss in the following our 
statement that weight functions (1 — s/M 2 ) k with large k are a bad choice for moments from 
the perturbation theory point of view, i.e. they are strongly nonperturbative. Indeed, the weight 
function in Eq. ([27]) has its maximum value at s max = M 2 l/{1 + k) (see Fig. ||). The integral 
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in Eq. (|28| ) is dominated by contributions from around this value s max for a smooth function 
p(s) as it is given by perturbation theory. The disadvantage of choosing such moments is that 
the factor (1 — s/M 2 ) k strongly enhances the infrared region of integration where perturbation 
theory is not valid. This can be seen by looking at the explicit convergence of the perturbation 
series for the moments (7|. 

To show this explicitly, we consider the integration of the running coupling with the modified 
weight functions in Eq. (|27|). The analogues of Eqs. (|21~1) , (f22j), and (|23| ) are necessary in order 
to find the behaviour of contributions due to logarithms (see Eq. (0)). The values of the 
coefficients of such contributions can be readily found in a concise form for arbitrary values of 
k at any giving finite order of perturbation theory (for a given power of the logarithm). For 
instance, the contribution of the term proportional to L = ln(M 2 /s) is given by 

where 7# is Euler's constant and ip{z) is the digamma function. In contrast to Eq. ([H]) the 
contributions in Eq. ( |25| ) increase as ln(fc) for large k. so the coefficients of the perturbation 
series for the moments grow relatively larger with increasing k. The contribution of the L 2 -term 
reads 

r M 2 ( a \ k / M 2 \ ^/e /fc+li\ 2 fc+1 i 

'I 1 -5?) to2 (^)lr(g7) +gH 7 * + ^+ 2))2+ ^ +2) (30) 

where ip'{z) is the first derivative of the digamma function. The coefficients in Eq. ([30]) grow as 
ln 2 (fc) for large k which has to be compared with Eq. (|22|). For higher powers of L the expressions 
become too lengthy to be presented here. The asymptotic behaviour of the perturbation theory 
coefficients for large k and fixed p is totally different from the behaviour obtained for the s l - 
moments, the direct moments from Eq. (fLOf) . 

The property of convergence of the resulting perturbation series for the moments can be 
reformulated in the language of effective scales for the moments themselves. This is similar 
to the approach where moments are reexpressed by moments Indeed, let us take the 

expansion of Eq. ( |T8| ) up to next-to-leading order, 

a M {s) = a M + PoLa 2 M + 0(a 3 M ) = a M + Poa 2 M ln(M 2 / s) + 0(a 3 M ) (31) 
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Figure 3: Different weight functions Wki(s) for (k,l) = (3,0), (2, 1), (1,2), and (0,3). 
Then for the s'-moments (i.e. the moments Mqi of Eq. (p8|) ) one has 

A) 2 



= % + 



l + l 



Ki + 0{a 3 M ) 



which translates into the effective scale M 2 e 1 '( l+1 > through the relation 



% + 



p a 



M 



l + l 



+ 0(4,) = a M (M r 2 e- 1 /('+ 1 )) + {a 3 M ). 



(32) 



(33) 



On the other hand, for the (M 2 — s) fe -moments given by Eq. (|28|) with large values for k and 
/ = the result of integrating Eq. ([31] ) reads 

m k0 = a M + A)<4 Mk) + 0(a 3 M ) = a M (M*/k) + 0(a 3 M ). (34) 

This leads to the effective scale M 2 /fc for large k. The results of these explicit calculations agree 
with a qualitative estimate based on the observation that the essential region of integration 
where integrals for the moments are saturated for reasonably smooth functions p(s) (which 
are just powers of logarithms in finite-order perturbation theory) is located around s max = 
M 2 Z/(Z + k). Obviously, the quantity a M (M 2 /fc) cannot be evaluated in perturbative QCD for 
large values of k. To be specific, for k = 2 we have 
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m 20 = a M + f3 a 2 M — + 0(a 3 M ) = a M 



and for k = 3 



6 



25 



(M T 2 e- n / 6 ) + 0( 



a M (0.16M 2 ) + 0(a 



(35) 



m 3Q = a M + A)4r^ + 0{a 3 M ) = a M (M T V 25/12 ) + 0{a 3 M ) 



a M (0.12M 2 ) + O(c4). (36) 
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Thus we find that ,in the language of effective scales, even a value of k — 3 requires the 
calculation of the effective coupling at a low scale which is definitely located in a nonperturbative 
regime. In this context we remind the reader that the effective scale Am for the effective coupling 
a M as defined in Eq. (|15|) is given by 

a m = exp(A; 1 //3o)A|s = 2.07A|g. (37) 



To fix the scale we take the value of A-^ = 350 MeV as determined from r decays P5| |. Then 
the process dependent parameter Am determining the effective coupling is numerically given 
by 

A 2 M = 2.07 x (0.35 GeV) 2 « 0.1 GeV 2 . (38) 

This is to be compared with the effective scale from Eqs. (|35| ) and (|3~6|). For M T = 1.777 GeV 
one finds 0.12 x M 2 = 0.36 GeV 2 which is definitely located in the nonperturbative region for 
the process under consideration. 

This is the situation with finite-order perturbation theory which is expected on the basis 
of asymptotic freedom of QCD and which is readily confirmed by a direct calculation. As the 
main output of this analysis we observe the movement of the related scale to where integrals are 
saturated and conclude that (M 2 — s) fc -moments for large k cannot be calculated in perturbation 
theory. Because the picture of explicit convergence due to the effects of running can be concealed 
by the use of a particular scheme (the MS-scheme, for instance) we have discussed the essence of 
the situation in an effective scheme where only effects of running are seen for the moments but 
no artificial constants related to the definition of an expansion parameter in a specific scheme 
occur. 

To conclude this section, we summarize our findings saying that demands from the theory 
side, basically dictated by the calculability in perturbation theory, and from the experiment 
side, basically dictated by the precision of experimental data, lead to conflicting requirements 
for the choice of weight functions as for where they should give the dominant contribution. 
Thus, one faces the usual conflict between experimental and theoretical accuracy which in 
our case is reflected by the range of (k, /)-values for the modified moments that are chosen as 
optimal observables. Having explicit perturbation theory formulae at hand (see e.g. Eq. (|T2|), 
(0)) (HP ; (HP ) (0)> an d (HID) one can establish the ultimate theoretical accuracy implied by 
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the asymptotic character of the perturbation series for a given experimental observable with 
any stated precision. This allows one to conclude which error - experimental or theoretical - 
dominates the uncertainty of an observable related to r decay physics. From the point of view 
of perturbation theory, large values of I are preferable while large values of k cannot be used. 

However, perturbation theory is not the end of the story. Instead, if the convergence is not 
satisfactory, this is an indicator that something is going wrong. There are contributions beyond 
perturbation theory, so if perturbation theory convergence is slow, perturbation theory itself 
is not reliable and nonperturbative terms become important. Nonperturbative terms for two- 
point correlators can be analyzed within OPE which can (and people believe, does) describe 
the spectrum at low energies better than ordinary perturbation theory. 

4 Power corrections 

In this section we consider the contribution of power corrections within OPE to the systems 
of moments given in Eqs. (|T0D and (p8|), and the interplay between the magnitude of this 
contribution and the structure of the perturbation theory series. In doing so we neglect, as 
usual, the weak ln(Q 2 ) dependence of the coefficient functions of local operators within OPE 
stemming from the anomalous dimension of the local operators to see the net effect of a power- 
type behaviour. Though this is a rather common practice in phenomenological applications, 
we do it for simplicity only; this dependence can be readily taken into account (see e.g. the full 
analysis of a realistic case below). The standard power corrections due to nonvanishing vacuum 
expectation values of local operators within OPE are relatively small for r decays but can play 
a crucial role for particular observables determined by special weight functions that lead to 
nonperturbative moments. Indeed, we saw in the previous section that (M 2 — s) fc -moments are 
not computable in perturbation theory for large values of k. Here we will see that one needs all 
power corrections which are known (or at least many terms) to bring the experimental results 
for the (M 2 — s) fc -moments in consistency with the theoretical framework based on OPE, i.e. 
the expansion in power corrections. 

For the system of moments in Eq. ( |10|) the contribution of power corrections reduces to a 
single term of the form (A 2 /M 2 ) 1 which decreases with / for A < M T where A is a typical scale 
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of the power corrections related to the nonperturbative scale of QCD. This restriction for A is 
a necessary condition for the applicability of perturbation theory. In fact, for A 2 /M 2 < 1 the 
perturbation theory contribution dominates in the total result for the direct s'-moments, being 
a necessary condition for using the technique. It is not quite clear, however, which numerical 
value should be chosen for A. The choice A ~ Aqcd is just too vague. In practice, a better 
choice for the scale could be A ~ m p or the mass of the respective low-lying resonance in the 
corresponding channel as for example the mass of the a\ meson for the axial correlator (see the 
discussion below). 

The perturbative contribution to the s z -moments with large values of I is saturated by high- 
energy contributions and therefore converges perturbatively, the convergence becomes even 
better with increasing I. Therefore, the s'-moments are theoretically computable with a strict 
control over the precision in perturbation theory. On the contrary, for I ~ and large values of 
k the system of mixed moments in Eq. (|28|) (which we call large k moments in the following) 
is saturated by low-energy contributions, i.e. basically by the contribution of the ground-state 
resonance, the existence of which is a common feature of the low-energy part of hadronic spectra. 
The large k moments are therefore completely nonperturbative. Power corrections to the large 
k moments come from many terms of OPE with local operators of high dimensionality. While 
these contributions are certainly decisive for reproducing the experimental spectrum, nothing 
definite can be said theoretically about such a sum of power corrections in any realistic case 
because high-dimensional power corrections are completely unknown numerically. 

An instructive example for the importance of power corrections to mixed moments with 
large k needed to reproduce the experimental spectrum of hadrons at low energies is the com- 
parison of moments for vector and axial channels. In the massless limit the perturbation theory 
series for the spectral density is the same both for the vector and axial channel leading to iden- 
tical expressions for the large k moments in perturbation theory. Experimentally, the large k 
moments are saturated by the contributions of the first resonances which are completely dif- 
ferent for both channels: the pion as a Goldstone boson in the axial channel and the p meson 
in the vector channel. Therefore, no method of summation of perturbation theory series alone 
can bring the perturbation theory results for the mixed moments with large k in agreement 



21 



with the experiment: if the method of summation is a regular one, the result is the same for 
both channels because it sums the same initial series. In this case perturbation theory for large 
k moments is in trouble and the power corrections have to provide the correct result. Because 
the asymptotic regime of very large k requires the consideration of all power corrections, this 
example shows that large k moments as in Eq. (^) cannot be used within the perturbation 
theory framework even if they are preferable from the experimental point of view. On the 
contrary, the system of direct moments in Eq. fllPf ) (also called I moments) can be calculated 
reliably in perturbation theory. These moments are saturated by perturbation theory which 
is also valid phenomenologically and known since long ago (chiral invariance, Weinberg sum 
rules) while the contribution of power corrections in this case is small and provides a necessary 
fine-tuning for better accuracy. 

Another way to understand what has been said up to now is to state that mixed moments 
with large values for k resolve the point-by-point structure of the spectral density at the origin 
which is definitely nonperturbative while mixed moments for large values for I resolve the point- 
by-point structure of the spectral density at large values for s where perturbation theory is more 
reliable and the standard asymptotic analysis is possible. 

After these preliminary remarks we give a realistic example where the magnitude of power 
corrections and perturbative terms can be seen quantitatively. In the consideration below we 
use the vector channel for the e + e~ annihilation into hadrons as an example which is relevant 
for the calculation of the electromagnetic coupling constant at Mz- This is the classical channel 
for vector mesons (p, uj). We give some details to elaborate the real situation in QCD within 
OPE. This example serves as a base for the next step - to simplify this picture to a model 
which contains only the gross features of the phenomenological structure of the spectrum but 
is simple and efficient for the purpose of investigating higher order power corrections within 
OPE. 

4.1 The vector channel contribution 

In this subsection we discuss the relative magnitude of perturbative and power corrections for 
a particular quantity related to the two-point correlator of hadronic currents of light quarks. 
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The quantity which is chosen is the value of the vacuum polarization function at the origin. 
This quantity is infrared (IR) sensitive, a fact that allows one to see the interplay between 
perturbative and power corrections. 

For light quarks the perturbative part of the vector (and axial) correlator is calculable for 
large values of Q 2 = —q 2 within the Euclidean domain. For the MS renormalization scheme it 
is given by 

n.^ V ,^) = ln(|)4 + ^(ln(|) + |-4f(3)) (39, 
/«s\V9 l2 /V\ /299 >/ ,\ ( a 2 \ 34525 715, . 25 >/K ,\ 

Eq. (|39|) is written for n/ = 3 active light quarks with the effective coupling a s = a^(fx 2 ). One 
is interested in an estimate for the value n llght (/i 2 , 0). However, II llght (/x 2 , Q 2 — > 0) cannot be 
obtained from Eq. (^) because there is no scale for light quarks, one can therefore not avail of 
a perturbative expression. Because singularities at small momenta are related to IR problems, 
it suffices to modify only the IR structure of the correlator IT llght (/z 2 , Q 2 ) in order to obtain an 
expression for all values of Q 2 . It is convenient to modify just the contribution of low-energy 
states to the correlator by using the dispersion relation 

n light (Q 2 ) = r (40) 

Jo s + Q 2 

where dimensional regularization is understood to be used for p llght (s). In fact, Eq. ( |40|) can be 
used for the bare quantities n llght (Q 2 ) and p llght (s). The modification should be done locally 
(i.e. with a finite support in the squared energy variable s in Eq. (0)) without changing the 
perturbative behaviour, in order not to affect any ultraviolet (UV) properties (/i 2 dependence) 
of the correlator Il llght (p 2 , Q 2 ). This requirement is important for retaining the renormalization 
group properties of the correlator. The low-energy modification of the perturbative expression 
for the spectrum is inspired by experiment: at low energies there is a well-pronounced bound 
state as a result of the strong interaction between quarks. Therefore, we adopt a model of 
the IR modification where the high-energy tail of the integral in Eq. fl40|) is computed within 
perturbation theory retaining the renormalization group structure of the result while in the 
low-energy domain there is a contribution of a single resonance. 
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For a generic light quark correlator in the massless perturbative approximation one intro- 
duces the IR modification 

P hsht (s) - p^ od (s) = F R 6(s - m R ) + p^(s)6(s - s ) (41) 

where Fr, tur, and sq are IR parameters of the spectrum. Note that they may not correspond 
to the actual experimental numbers. Substituting the IR modified spectrum in Eq. ([4lD into 
Eq. (0), one finds 

n?Il> 2 , 0) = 4 + "> f £) + \ + " fin f f) + g " 4C(3)) (42) 
m R \s J 3 7T \ \ s J 12 y 

We identify rriR with a mass of the low-lying resonance which is the only input giving a scale 
to the problem. The IR modifying parameters Fr and So are fixed from quark-hadron duality 
arguments. Notice the 0(a^) difference between Eq. (|39|) and Eq. (^) given by the additional 
term — 37i~ 2 /8. This is the so-called "7r 2 -correction" (see e.g. Ref. |39[]) which can be rewritten 
through C(2) = tt 2 /6. 

Next we use the OPE with power corrections that semi-phenomenologically encode the 
information about the low-energy domain of the spectrum through the vacuum condensates of 



local gauge invariant operators |4(| in order to supplement the perturbative contribution. The 



OPE for the light quark correlator reads 



n& 2 , Q 2 ) = n£ od (/A q 2 ) + <^> + <p> + . . . (43) 



The quantities (04,6) denote nonperturbative contributions of dimension-four and dimension-six 
vacuum condensates. These contributions are UV soft (they do not change the short distance 
properties of the correlator) and are related to the IR modification of the spectrum. For the 
purposes of fixing the numerical values of the parameters Fr and so which describe the IR 
modification of the spectrum one needs only the first two power corrections 1/Q 2 and 1/Q 4 . 
Note, however, that the coefficient of the 1/Q 2 correction vanishes because there are no gauge 
invariant dimension- two operators in the massless limit. Computing the IR modified and OPE 
supplemented correlation function, we find finite energy sum rules (FESR) that alow one to fix 
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the parameters Fr and s JIT] 



Fr = So |i + ^+^)V^1b^ +^ + ^1 + 0(02), 

^ = Iji+v+v^^^+^ + f)}-^ 4 ^ ^- (44) 

We treat (O4) as a small correction and take its coefficient function as a constant (the total 
contribution is renormalization group invariant). Eqs. ( pH| ) fix Fr and Sq through m R and 
(O4). In using higher order terms in the OPE expansion (for instance, (Oq)/Q 6 ) one could 
avoid substituting m 2 R from experiment because within the IR modification given in Eq. ( [41| ) 



the IR scale is determined by the dimension-six vacuum condensate (Oq) [41|. We do not do this 
because the primary purpose of the present analysis is to find the low-scale normalization for 
the electromagnetic coupling and not to describe the spectrum in the low-energy domain. The 
use of the experimental value for the resonance mass m R makes the calculation more precise 
because the numerical value for the (Oq) condensate is not known well (cf. Ref. |[42|| ) . 

The leading order solution for Eqs. ( [44]) (upon neglecting the perturbative and nonpertur- 
bative corrections) is given by the parton model result Sq = 2m 2 R and Fr = Sq = 2m R which 
is rather precise. This solution has been used for predicting masses and residues of the radial 
excitations of vector mesons within the local duality approach where the experimental spec- 
trum is approximated by a sequence of infinitely narrow resonances ||43|| . The condensate of 
dimension-four operators for light quarks is given by 

(O,) = ^(l+ 7 -^) (-G 2 ) + 2tt 2 ( 1 + (m u + m d ){(uu) + (dd)). (45) 

3 V D7T / 7T \ 671/ 



The order a s corrections to the gluon condensate were calculated in Ref. H44| . We retain small 



corrections proportional to the light quark masses and treat them according to the approxima- 
tion of isotopic symmetry for the light quark condensates, (uu) = (dd), which is rather precise 
for u and d quarks. The quark condensate part of Eq. ( f45| ) is given by the partially conserved 
axial current (PCAC) relation for the n meson, 

(m u + m d ) (uu + dd) = -flm\. (46) 

Here f n = 133 MeV is the decay constant of the charged pion and m n = 139.6 MeV is the mass 
of the charged pion. For the numerical standard value ((a s /n)G 2 ) = 0.012 GeV 4 of the gluon 
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Figure 4: Spectral densities for the models for the vector channel (corresponding to Eq. 
and the axial- vector channel (corresponding to Eq. (|57|). The arrows indicate the narrow width 
resonances, given by distributions. 



condensate j|0| and a s /ir = 0.1 one finds 



(O a ) = \ (l + ^) <^G 2 > - 2tt 2 (l + ^) flml = 0.037 GeV 4 . (47) 

The relation sq = 2m 2 , where m p = 768.5 MeV is a mass of the lowest (p meson) resonance 
in the non-strange isotopic 1 = 1 vector channel, is rather precise numerically. The gluon 
condensate gives a small correction to the basic duality relation s = 2m|, for light quarks. 
The same is true for the 1 = channel where the lowest resonance is the u meson with a mass 
m tu — 781.94 MeV. So we conclude that the model in Eq. ( f4~T|) is rather accurate if the scale is 
fixed by the low-lying resonance. 



4.2 Models for the vector and axial-vector channel 

Based on the above example we investigate models for non-perturbative effects resulting from 
standard power corrections [0 that allow us to obtain all power corrections in a concise form. 
We adhere to models for the vector and axial channels which are simple enough to be calculable 
and which retain the main features of the spectra. These models are depicted in Fig. [|. With 
these models at hand we have all power corrections available for the discussion of their role in 
moments for the spectral density in different channels. 
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As a model for the vector channel we take the spectral density of the previous subsection 
in a simplified form (m = my), 

p v (s) = 2m 2 5(s -m 2 ) + 6(s -2m 2 ). (48) 

By using the dispersion relation we obtain the correlation function 

2m 2 ( u 2 \ 

n "« 2) = + ln J + (49) 

with necessary subtractions. The expression in Eq. Q4"5|) is used to generate all power corrections. 
Indeed, for Q 2 ^> m 2 one finds 

n^)^) + £(-fH^) 

where the first term is the leading order perturbative contribution in a s and the remaining 
terms are power corrections. Note that the analytic properties of the expansion in Eq. flSTp up 
to any finite order are different from the exact (though model-dependent) result in Eq. fl4~9|). 
This is rather a common feature: analytic properties of approximations for the correlators can 
be different from those of the exact result. In some instances this restricts the precision and 
may lead to a misuse of approximations in areas where they do not work. 
Writing the OPE for the correlator in the general form 

nv(Q 2 ) = m(g + £^p <«) 

for the model in Eq. ( f49|) we find 

-2m 2 ) n . (52) 



The first two coefficients vanish, c\ = C2 = 0. The vanishing of c\ is in full agreement with 
the fact that there exist no dimension-two operators in realistic cases while the vanishing of c-i 
means that the gluon condensate is neglected in this model (which is justified numerically for 
the realistic case of r decays, see also the previous section). Note that these two constraints are 
built-in requirements for our models. We just constructed the models in this way. The third 
coefficient reads 

2 

c 3 = -g 7 ™ 6 - (53) 
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Phenomenologically this coefficient is related to the vacuum expectation value of local four- 
quark operators which in factorized approximation is given by 

<<#> = - 8 ~^a s m 2 . (54) 



\1 

Its numerical value is approximated reasonably well by the expression in Eq. (|53| ) with m = m p . 

For further comparison we give the value of the dimension-six term in the axial channel, 

imA , 14087T 3 , l9 , . 

(0£) = -^-a s (qq) 2 . (55) 

This vacuum expectation value is larger and has a sign opposite to that of the vector channel. 
The first few terms of the model OPE in the vector channel read explicitely 



> 2 \ 2m 6 2m 8 22m 10 26m 12 

g 2 " ) ~ W + ~W ~ "5Q 1 " + ~W 2 



^v(Q 2 )=ln ^ -— K + ^-—^ w + - 7 ^ + ... (56) 



The term of dimension eight is still available in OPE as expressed through the vacuum expecta- 
tion value of local operators, even though it is very poorly known numerically |45j, while higher 



order terms were never used in phenomenological applications. The expression in Eq. (|50f) 
shows also the actual scale of the expansion in the vector channel, Sq = 2m 2 . The scale A 
with A ~ A QCD ~ A^jg ~ 350 + 400 MeV or A 2 ~ 0.25m 2 does not fit the scale of the power 
corrections in this model. 

There is a numerical cancellation between the resonance and continuum contributions to the 
coefficients c n of the first several terms while for higher order terms (large values of n) the scale 
So = 2m 2 dominates. This cancellation is one of the reasons for the success of the Borel sum 
rules for the p meson in the vector channel. One can go very low in the Borel parameter M and 
still finds power corrections to be small if only power corrections up to dimension-six operators 
are included (with a nonvanishing gluon condensate there is also an additional cancellation 
between the dimension- four and dimension-six terms). Note that the expansion for Adler's 
function D(Q 2 ) (or for the renormalized correlator function U(Q 2 )) converges for Q 2 > 2m 2 . 
This is related to the simplicity of the analytic structure of the correlator in the model and 
therefore to the simplicity of the model spectrum. Therefore, if one is to believe that the shape 
of the spectrum is really close to that of the model in Eq. fl4"8]) , the scale of power corrections 
should be close to the value s = 2m 2 . We discuss this issue later on in more detail for a variety 
of different models. 
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Next we consider the axial part of the correlator or the axial channel. Because of the 
presence of the pion, in this case the spectrum at low energies is drastically different from the 
one for the vector channel. All axial-vector resonances (with spin 1) have a finite mass. In the 
massless limit there is theoretically a Goldstone mode - corresponding to the observed pion - 
with spin zero contributing to the correlator of the axial-vector current (this is the reason why 
we call it more generally axial correlator or axial channel). The main mass scale is the mass 
m ai of the axial-vector meson a\ which we express by rr? a = 2m\ for further convenience. The 
underlying reason for this is Weinberg's relation: at some point one can identify the scale 
with the scale in the vector channel which is given by the p meson mass m p . A possible model 
for the spectrum in the axial channel reads 

p A (s) = m 2 A 5(s) + m A 8(s - 2m 2 A ) + 6(s - 2m 2 A ) (57) 

where the first term is the pion contribution, the second one is contribution of the a>i meson, 
and the third represents the continuum. There is no gap left between the second resonance and 
the continuum. The correlator in the axial channel is given by 

the expansion at large Q 2 reads 



n A (Q 2 ) = 




where 8 n \ is the Kronecker symbol. Here the contribution of the dimension-four operator is 
again zero while the dimension-six contribution is positive and larger than that in the vector 
channel, which is the case also in the (model independent) OPE (see Eq. (|55|) ). While the 
continuum contribution (the logarithm and the 1/n part in the sum in Eq. (|50|)) remains the 
same, the factor — l/2 n_1 in case of the vector channel is replaced by —1/2 (for n > 1) in 
Eq. (|59|). Therefore, higher order power corrections for the vector channel are dominated by 
the continuum while for the axial channel they are dominated by the resonance contributions 
and are generally larger - the mass of the a\ meson gives the scale both for the resonance 
contributions and the continuum threshold in this particular model. However, this cannot be 
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quantitatively checked at present because the numerical values of the higher order condensates 
are not known phenomenologically with sufficient accuracy. 

It is instructive to compare the results for the model OPE expansions in the two channels, 

TT <f*\ In 2m " V 4- 22m V + 26m V J 

MQ) - ln l v giJ-^ + ^r-^gir + ^ + --- 

„ , f u 2 \ km\ Am\ 48m 1 , 64m 1 ? , . 

Ka(Q 2 ) = fn ^ )+-A-—A + --^--—^ + ... (60) 



The coefficients of the power corrections are different and the signs of corresponding terms are 
opposite. A numerical comparison is possible if one identifies the scales m^ = my = m in both 
channels which is a rather reasonable phenomenological approximation because of Weinberg's 
relation m? ai = 2m 2 p . However, one should keep in mind that both models only capture the 
gross features of the spectra while the fine details (visible at high resolution in the energy) 
are different and of importance for the numerical magnitude of higher order power corrections. 
Indeed, higher order power corrections are sensitive to the fine details of the spectra and 
cannot be reliably determined within rough models. Therefore, while the models proposed in 
this subsection allow one to calculate power corrections up to any order, a reasonable accuracy 
is expected only for low orders in 1/Q 2 - high-order power corrections can resolve the fine 
structure of the spectrum which is not caught by the gross models. It is expected that for the 
first few terms of the power expansion the accuracy is rather good while for high-order terms 
it can be only an order of magnitude approximation. Still such simple models based on the 
gross features of the spectrum are definitely useful for a general analysis. Higher order terms 
are more sensitive to details of the spectrum and can be predicted only with large errors. We 
discuss these issues below. 

The moments corresponding to the model spectral density for the vector channel for k = 
read (assuming 2m 2 < M 2 ) 



(2m 2 \ l+X ( l + 

One observes that the perturbative contribution is represented by the first term on the right 
hand side of Eq. (|6lD. The power corrections are given in what follows. The combined per- 
turbative and power correction structure is a natural order for the direct s'-moments. For 
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large I the contribution of the power corrections decreases and the moments are saturated by 
perturbation theory, i.e. if m 2 <C M 2 , the power corrections for the moments Mqi die out fast. 
For mixed moments with I = and arbitrary k we obtain (again for 2m 2 < M 2 ) 

M'-sr+^SM)'- (62) 

The magnitude of these moments tends to zero for large values of k and definitely cannot be 
represented perturbatively. Indeed, according to the structure of OPE one should represent the 
moments as 

M k0 = 1 + A fc0 (63) 

where A^o gives the contribution of the power correction terms which are usually considered to 
be small. For large values of k, however, one has M k o — » (as Eq. (B3) shows) and A^o — » —1 
from Eq. (|B3|), i.e. the moments are not given by an expansion near the perturbation theory 
result, M fc0 = 1. The decomposition in Eq. ( |63"D is therefore useless for large values of k because 
Afco is not small. For the first few orders in k one has 



/ TV? \ ( TV? \ ( TV? 

A 00 = Aio = 0, A 20 = -2 ( — ) , A 30 = -8 | — ) + 8 | — | 



M 2 / ' \M 2 ) \Ml ) ' 

/m 2 \ 3 /m 2 \ 4 /m 2 \ 5 

A « = - 20 fe) +40 (m?) - 22 U) ' 

( m 2 \ ( m 2 \ ( m 2 \ ( m 2 \ 

A » = - A0 {m) +120 (m?) - 132 fe) +52 (m;J 

Aw . _ 7 og) 3 + 28 og) 4 - 4 a 2 (|) 5 + 3 6 4 g) 6 - 114 (g) 7 . , M) 

If m 2 -C M 2 , the power corrections for the moments M^o are still basically given by the lowest 
order term. If 2m 2 is close to M 2 as it is actually the case in r decays (the conclusion is based 
on the model spectrum), power corrections for the moments M i are given by a single operator 
of dimension / and are relatively small while for the moments the power corrections are 
given by a linear combination of all operators up to the specified order. It is instructive to 
rewrite the k moments in terms of the actual scale 2m 2 . One has 

a a a 1 /2m 2 \ 3 /2m 2 \ 3 1 /2m 2 \ 4 

A °° - Al ° - °' A2 ° - 'lyM 2 ) ' Aso " AW J + l\Mf) ' 

5 /2m 2 \ 3 5 /2m 2 \ 4 11 /2m 2 \ 5 
A 4 o = -7T Ktt + 



2 V Ml } 2 V M 2 J 16 V M 2 
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/2m 2 \ 3 15 (2nf_Y _ 33 (2rrf\ 5 13 /2m 2 \ 6 

35 /2m 2 \ 3 35 /2m 2 \ 4 231 /2m 2 \ 5 91 /2m 2 \ 6 57 /2m 2 \ 7 
Ae0 " ~T\M?) + Y\M?) ~1g[m*) + 16\-M?) ~ 64 \~M~f J ' (65) 

Looking at Eq. (|65| ) one can conclude that keeping only the first contribution can give a com- 
pletely wrong answer for the total result for A^o at large k. Therefore, being expressed through 
the natural scale, i.e. the scale where the analytic properties of the correlator change, the two 
sets of moments behave differently. The direct moments have power corrections in accordance 
with the expectation that these corrections are powers of the relevant scale while the numerical 
magnitude of mixed moments is non-predictable on general grounds without a detailed knowl- 
edge of the spectrum at low energies. Note again that we neglect the effects of the anomalous 
dimensions of local operators in our simplified consideration based on a model spectrum which, 
however, is also a common practice in phenomenological applications. 

The dominance of the resonances for the axial channel becomes obvious if we consider the 
moments. For 2m\ < M 2 the direct moments are given by 

M 00 = l, Moz = l+^#V + 7l-i±iV (66) 



Ml 



while the mixed moments for Z = read 



M M 4-fp<* +1) g +(fc+ 4-fyf. (67 , 

Because of the second term in Mko which comes from the pion resonance at s = 0, the moments 
Mfco will not vanish for increasing values of k but will increase linearly in k. This is a feature 
worth emphasizing. For higher and higher values of k the moments M^o in Eq. (|62[) go to zero 
if used for smooth spectral functions on the cut which is the case if the spectral densities are 
calculated in perturbation theory. However, the experimental shape of the spectrum is different, 
and the spectral density partially behaves like a distribution. The total pion contribution for 
example is concentrated at one point which is not affected by the decrease of the essential 
support of the measure with increasing k. Therefore, due to the normalization chosen the pion 
contribution grows with k. These two features of the experimental k moments, namely the 
decrease of the vector contribution and the increase of the axial contribution due to the pion 
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Table 1: Moments Mu for k = for the vector cannel and the axial channel model, using the 
exact result M^f as well as the power series expansion up to the third order term in 1/Q 2 , Mffi. 
We use the estimate m 2 A = my = m 2 p = 769.3 MeV @. 



(as found in Eqs. (62) and (pf\l) are very essential for the comparison with the theoretical 
description of the r-moments to be successful. Indeed, one has 

(the normalization factor Nfa is chosen so that Mfc = 1 for p(s) = 1). Remark that the 
normalized r lepton decay rate R T as in Eq. ([/]) is just Mq . The theoretical expression for the 
moments changes slowly with k, therefore the sum of vector and axial parts behave better than 
these two parts separately. 

With the explicit models in Eqs. ( [48|) and fl57"D at hand we can make a quantitative analysis. 
We assume that these models give the exact results at any order of the power expansion. Then 
we consider the moments and figure out when it is enough to keep only the three terms usually 
available in phenomenology in order to have a reasonable (given) accuracy. Table [l] shows the 
situation for the I moments. One sees that the perturbation theory contribution dominates the 
results. For both vector and axial correlators the accuracy of the three-term approximation is 
better than 2% and improves for large values of I as expected. 

Table |2| shows the corresponding situation for the k moments. Here one sees the dominance 
of nonperturbative contributions to the vector and axial correlator results. For k = 4 the 
accuracy is already about 10% and deteriorates fast. For the sum of vector and axial correlators 
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Table 2: Moments M^i for I = for the vector cannel and the axial channel model, using the 
exact result as well as the power series expansion up to the third order term in 1/Q 2 , . 
We use the estimate m\ = my = m 2 p = 769.3 MeV 0. 
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Table 3: Relative deviation of the moments M k i for the power series up to the third order 
from the moments M kl for the full result containing the vector and the axial-vector channel 
model. We use the estimate m 2 A = m 2 B = m 2 = 769.3 MeV 0. 
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the accuracy of the three-term approximation is better as one can see from the last line of 
Table [| This proves our statement about the behaviour of the sum of vector and axial part 
which can also be seen from the expansions in Eqs. fl60|). It is known up to the third order 
term from phenomenology where the corresponding values are given in Eqs. ([54]) and (0). Note, 
however, that the phenomenological results are obtained only in the factorization approximation 
for the vacuum expectation values of the four-quark operators. Thus, the pion contribution 
is crucial in keeping a reasonable accuracy for the sum of axial and vector correlators for the 
large k moments. This corresponds to the local duality picture. 

Our general arguments are becoming more transparent if we consider mixed moments M^i 
with non-zero values of /. This is done in Table [5] where we give ratios 

for the r-moments defined in Eq. (|68|). Then only r-moments for large values of I and small 
values of k show perturbative behaviour. The bottom line in Table |3| is drastically different 
from the rest of the table, the reason being that the massless pion only contributes to moments 
with / = 0. This makes the large k moments still reasonably precise using only third-order 
power corrections. 

4.3 Modifications of the model for the axial channel 

Here we discuss the choice of the spectrum in the axial channel. One can suggest forms differing 
from those given in Eq. flSTp. The reason for this freedom is that we have two resonances which 
are essential at low energies - the pion and ai meson - that bring in four free parameters: the 
two residues of the resonances, the mass of the a\ meson, and the beginning of the continuum. 
Therefore, two requirements for the vanishing of the power corrections 1/Q 2 and l/Q 4 are not 
sufficient to uniquely fix the low-energy spectrum in a simple way (with only one scale). There 
can be many different additional requirements all reasonably close to the experimental data. 
Note that the a\ meson shows up as a broad resonance which provides additional freedom in 
parameterizing its contribution to the spectrum. For a precision application the best choice is 
to fix the pion residue and the mass of the a\ meson to their experimental values while leaving 
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Figure 5: Spectral density for the modified axial- vector channel model given by Eq. (|72|) . 

its residue and the continuum threshold to be determined from OPE. In fact, the experimental 
value of /„- and simple duality requirements give a rather good description of the axial channel. 

Another possibility is to fix the residue of the a\ meson to be 2m 2 A as in the vector channel 
(see Eq. (^8|)). Then the spectrum in the axial- vector channel reads 



Pa 



» = 2( v / 2 - l)m 2 A 8{s) + 2m 2 A 5{s - 2m 2 A ) + 6(s - 2 v / 2m^) (70) 



which also leads to a result in which the coefficients of the first two power corrections vanish. 
The third coefficient of the large Q 2 expansion of the axial correlator now reads 

8 2 
c^ = -(3-2^) = 0.457528 < - (71) 
3 3 

(cf. Eq. (^0])) which is smaller than in the model given by Eq. (|5"TD. It is also smaller than in 
the vector channel which can lead to a (slight) contradiction with the phenomenology if the 
scales my and uia are identified and factorization of the four-quark condensate is taken as 
being exact. Still it is numerically acceptable from the phenomeno logical point of view. 

Having fixed the mass of the a>i meson to 2m 2 A one has a one-parameter family of spectra 
in the axial-vector channel 

PA (s) = An 2 f 2 5(s) + An 2 f 2 a 5(s - 2m 2 A ) + 6(s - s ) (72) 

with sq being a free parameter. Then 



The spectral density of is modified model for the axial channel is shown in Fig. |5|. One constraint 
for the construction of the model is Sq > 2m A . However, this constraint is quite weak. One 
can actually admit a resonance as a bump in the continuum which is almost the case for 
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the ax meson. The positivity constraint for f 2 is given by sq < Am\ and has to be taken 
seriously. Note that at so = Am\ the pion decouples and does not contribute (f % = 0) to 
the spectrum of the correlator which is not the case in reality. For some values of So in the 
interval 2m\ < sq < 4m 2 A the spectrum can be fixed by taking the experimental value for 
/„■. A possible additional constraint is given by the value of the dimension-six condensate. 
This was used for the finite-energy sum rule analysis and turned out to be rather successful 
phenomenologically [f|I|. 



The residue of the a\ meson and the continuum threshold for the model in Eq. ( ff^) (as 
being expressed in terms of the pion decay constant /„. and the mass m ai of the a\ meson) are 
given by 



s = 2m\ (l + y/l - 4n*f*/m? A ) , An 2 f 2 a = s - An 2 fl (74) 
The consistency constraint reads 

^ 2 f 2 Jm\ < 1. (75) 

Experimentally we have ^ 2 f 2 ~ m 2 A which is close to the model without a gap. Our conclusions 
about moments are valid for these models as well. 

The considerations obtained for the simple models lead to the conclusion that the structure 
of power corrections for the moments with weight (1 — s/M 2 ) k is strongly infrared sensitive 
and nonperturbative. One can expect that this conclusion will be retained for a realistic model 
(QCD), i.e. that power corrections are very important and that the perturbation theory con- 
tribution is strongly suppressed for the moments with weight (1 — s/M 2 ) k . This was to be 
expected because these moments are saturated by resonances. 



4.4 A note on theoretical models for power corrections 

For the theoretical analysis of large k moments one needs all power corrections. These are 
available only in specific models. In the previous parts of this section we have discussed models 
for the spectra based on the experimentally measured shape of the spectra while taking into 
account only main characteristic features of them. However, there are also theoretical results 
for the spectra that are usually obtained in simplified models of quantum field theory and 
which, as one believes, can have some relation to the low-energy hadron spectra in QCD. 
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Figure 6: Spectral density of the narrow resonance models for the vector and axial-vector 
channel as given in Eqs. ( f78|) and (|82|). 



The old example is two-dimensional QED (QED2) which was solved by Schwinger f46fl . The 
OPE for current correlators in QED 2 was analyzed in Ref. |47| . In this model the scale of 
power corrections is given by the dimensional charge e (the model is super-renormalizable). 
Both vector and axial- vector currents have the same scale e 2 /7r, in fact they are given by a free 
massive field and are trivial in this respect. The photon acquires a mass due to the anomaly. 
The interesting channels are that of the scalar and pseudoscalar currents, the correlators for 



which are known exactly and give an example of a really nontrivial spectrum [47 



There are some results concerning all-order perturbation theory diagrams calculations in a 



ladder approximation | 4q| . These results can be used for constructing realistic four-dimensional 
models where all power corrections are explicitly known. The sum rule analysis of OPE in such 
models can reveal some features different from the standard phenomenlogy [49|. 

A popular model is the 1/N C approximation both in four and two-dimensional Abelian and 
non-Abelian theories (see e.g. Ref. p^, |51|) where the spectrum is given by infinitely narrow 
resonances, as shown in Fig. |6|. This type of spectrum is predicted by a Regge analysis or by 
local duality arguments. A recent analysis along this line was presented in Refs. |52| , |53| . 

The model of narrow resonances inspired by the 't Hooft model was analyzed in QCD on 
the basis of local duality [^3|. The spectrum was studied within the local duality approach 
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where one has |H3 



m 2 = (2n + l)m 2 , f 2 = 2m 2 n = 0, 1, . . . (76) 

for the vector channel and where m p is the mass of the the ground state p meson. For the 
axial-vector channel including a massless pion as Goldstone boson (axial channel) the result of 
the local duality approach coupled with the 1/N C approximation reads j|3| 

1 
2 



m 2 n = nm 2 ai , f 2 = 4vr 2 / 2 = X = m% f 2 = 2m 2 p = m\ , n = 1, 2 . . . (77) 



where m ai is the mass of the ground state a\ meson. The structure of the spectrum reflects the 
classical results on chiral symmetry and Weinberg's relations for axial-vector and vector meson 
masses which is realized if one identifies the scales in both channels, m 2 = 2m 2 = 2m 2 . This 

; til P 

identification leads to a simplified picture where two chains of resonances are simply shifted by 
an amount m 2 = m 2 . This is the gross structure of the spectrum. Experimental results differ 



quite a bit from this picture, but such details can be accounted for by using the OPE glj . 

Within the model based on the local duality approach the summation of all resonances 
results in the vector correlator 

oo f2 oo Q 2 

n.(Q 2 ) = E^tts = E 



n=0 ml + Q 2 n f (2n + l)m 2 + Q 2 

00 1 (Q 2 + m 2 \ 

= — — r— 77 = — ib — + subtractions. (78) 

^o n + (Q + m 2 )/(2m 2 ) y \ 2m 2 J K ' 

Here ib(z) is the digamma function. To find the link to the OPE we consider the large Q 2 

behaviour of the model correlator in Eq. (|78|). The large z expansion of the digamma function 

ib(z) reads 

1 B 2k _ / 1 



*(-) = -»- s-gs^ + o^J P-) 

where S2A; are the Bernoulli numbers 

^= ( "ff' C(a) ' B2fe+1 = ° (80) 

and C(2^) i s the Riemann £ function. This expansion is an asymptotic expansion. For the 
correlator one then obtains 



. ., / u 2 \ m 4 m 6 7m 8 m 10 
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where the renormalization scale \i comes in through the subtraction term. Note that the 
expression in Eq. ( |3~1~D is quite different from the result obtained in Eq. (|50 ) . 



In the axial channel the model based on local duality (with parameters from Eq. (|77|)) leads 
to the expression 

oo p2 x2 oo ^2 



™l + Q 2 Q 2 ^ m 2 ai n + Q 2 
ml ^ 1 m\ ( Q 2 \ 

= W + S ^ov< = ~w - * l< J + subtractions (82) 

and the large Q 2 behaviour reads 

/Ox / u 2 \ mi mi . , 

n ^ 2 ) = ln (^) + i# + #-iw + # + - (83) 

In this case all power corrections of dimension 2(2k + 1) vanish because of I?2fc+i = 0. This 
contradicts the results of OPE since there is no explicit (symmetry) reason for such a property. 
At moderate orders of n the structure of the expansion of Eq. (|52D as given in Eq. ([83]) is incon- 
sistent with the asymptotic expansion expected from the OPE (see Eqs. ([48]), ([49]), (|50|), (|58|), 
and (|59|)) and, therefore, is not supported by phenomenology. Though the model of infinitely 
narrow resonances is rather attractive, this inconsistency was the reason that this model was 
not analyzed quantitatively in the literature and was instead considered only as a rough ap- 
proximation for the first few moments (duality relation). Indeed, the zeroth moment is correct 
(no 1/Q 2 terms) and the first moment is correct in the approximation of a vanishing gluon con- 
densate contribution but the vanishing of the dimension-six terms are definitely unacceptable 
phenomenologically. The high order 1/Q 2 behaviour can be easily corrected by shifting the 
parameters of the first one or two resonances. But then the simplicity of the functional form 
of this model with its one-parameter dependence for the whole spectrum is lost. 

The qualitative difference of the models given in Eqs. (|8T]) and (^) from the previous case 
with a continuum contribution, however, is the analytic structure of the correlators. In the 
narrow resonance model one only has a single dimensional parameter m 2 and one would expect 
the power corrections to behave as (m 2 /Q 2 ) n with m 2 determining the scale. However, the 
coefficients of the power corrections given by the Bernoulli numbers grow quite fast as the large 
k behaviour demonstrates, 

/ , \ 2fc 

B 2k = A(-l) k - 1 V^i— J (l + o(l)) (84) 
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(the notation o(x p ) indicates contributions which fall off to zero faster than x p for x — *■ 0). 
Thus, the expansion in Eqs. ( pi]) and fl83|) is only asymptotic. The reason for the divergent 
nature of the series of power corrections is that for the models given in Eqs. ( [Tip and ( |82] ) 
there are poles located arbitrarily far away from the origin Q 2 = 0. In this model, therefore, 
there is in fact an infinite number of scales nm? for positive integers n. This is the reason 
why perturbation theory does not work even at sufficiently large s. But if one approximates a 
chain of resonances by a continuum starting from some threshold Sq, one instead gets a finite 
radius of convergence of the order s , Q 2 > s and a well working perturbation theory picture 
at s > s . 

Let us compare the behaviour of moments for such models. Because of the complicated 
argument of the ip function in Eq. ( |78|) for the vector correlator, power corrections in the form 
of an expansion in 1/Q 2 are not directly given by the expansion of the ip function. This makes 
the calculation of an arbitrary term for the vector correlator clumsy. In case of the axial 
correlator this expansion is straightforward as one can see from Eq. fl82"|) . Therefore, we shall 
only analyze the axial channel but will give a functional relation that holds for both correlators 
in these models. Identifying 2m 2 with m 2 ai , m 2 ai = 2m 2 one finds 

2 

Tfl 

Tl A (Q 2 ) = Tl v (Q 2 -m 2 )-— (85) 
or, when rewritten in the Minkowskian domain, 

2 

Tfi 

U A (q 2 ) = U v (q 2 + m 2 ) + — . (86) 

This relation is inspired by particular models for the spectra but can be checked for experimental 
data as well (for Adler's functions which are independent of subtractions). 

The narrow resonance model for the vector and axial channel does not display the difference 
between low and high energies. The only criterion of perturbation theory calculability is the 
length of the averaging interval while its position in terms of energies is almost unimportant. 
This is a natural feature of the translation invariance of the spectra in these models. We now 
discuss the moments in the axial channel in more detail. In order to obtain simple estimates 
we take the averaging interval to be 5m 2 (in "real life" the averaging interval is given by M 2 = 
5.15m 2 but 5m 2 = 5m 2 is a good approximation). We check the quality of the decomposition 
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Table 4: Deviations A\ L and Ajy (see Eq. (0)) of the moments for the narrow resonance vector 
and axial channel model with respect to the parton contribution A\[ A = 1 for k = or / = 0. 



from the perturbation theory point of view, i.e. we write the result in the form 

M u = 1 + A kl . (87) 
The direct / moments obtained for the power expansion in Eq. (^) are (for / > 0) given by 



mi x 



M ol = l + (l + 1)1^1 B 2l . (88) 



For large values of / the decomposition as given in Eq. (|87D is useless because the coefficients 
are growing very fast and thus the second term becomes large (see Eq. @)). Still, for odd 
moments the decomposition is precise because then only one power correction is picked up and 
its contribution vanishes. However, the result based on the asymptotic value Mqi = 1 is quite 
different from the exact answer 

(S) 2,+ (f) 2 ' (89 » 

The reason for this is that the asymptotic behaviour is too crude to resolve the position of the 
poles near the origin. This can be seen from the expression for the vector channel where the 
appropiate expansion parameter for power corrections is 1/(Q 2 + m 2 ) instead of 1/Q 2 - In other 
words, the scale M 2 is two small to be able to reflect any information about remote poles which 
give contributions to the asymptotic behaviour at large Q 2 . 

For the mixed moments in the axial channel our result are given in the lower part of 
Table |j. We give only numerical results because it is not straightforward to represent the 
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decomposition analytically. The message is clear: both systems of moments - direct and mixed 
- are nonperturbative which is expected for this model where for any energy region the spectrum 
has nonperturbative character (close to a singularity). 

An even more symmetric picture emerges in the vector channel. In this case one does not 
see such a drastic difference between moments as should be expected from the perturbation 
theory point of view: the situation is symmetric, the scale of "bumpiness" is the same at the 
both ends of the spectra. This fact makes both systems of moments almost equivalent. If we 
take the average interval 6m 2 , the result for exact moments is the same for both systems of 
moments because of the relation 

(1 - 1/6)" + (1 - 3/6)" + (1 - 5/6)" = (1/6)" + (3/6)" + (5/6)", (90) 

as shown in the upper part of Table f|. In general, for the vector current there is just an exact 
symmetry of the model spectrum s — *> (2km 2 — s) which makes the result of the calculation for 
both systems of moments simply equal. Still the contribution of power corrections is different 
for different systems of moments and, moreover, does not allow to restore the shape of the exact 
spectrum. 

This symmetry of the simplified model spectrum is definitely violated in the realistic phe- 
nomenological spectrum, making the region of high energy essentially different from the low 
energy domain. This is because of the asymptotic freedom in QCD where at large energies the 
spectrum becomes a more and more continuous function because the features of the resonances 
as distributions becomes less and less pronounced. This also means that the model for the 
spectrum as an infinite chain of infinitely narrow resonances does not properly incorporate the 
asymptotic freedom of QCD in a sense of scale invariance violation. 

We should note here that if there is a new situation at large energies (as the production of 
heavy quarks near their threshold) one should do the analysis as if the energy were counted 
from the respective thresholds. It is not the absolute value of the energy that matters but the 
distance from the strong interaction region. 

As a last remark we stress that even the formal knowledge of the asymptotic series in 
l/Q 2 does not allow one to restore the spectrum point-wise. This reflects the same situation 
as in phenomenological analyses, namely that the determination of the spectral density is 
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an improperly posed problem (inverse problem). In our particular example one can actually 
observe this feature. For the axial channel one can still figure out the spectrum from the 
coefficients by restoring the full function but for the vector channel it is not straightforward. 
There is actually no systematic method for doing so. Taking just the "leading" asymptotics for 
the coefficients one can easily obtain identical results for the axial and vector channels which is 
definitely an unwelcome result as regards low-energy phenomenology. Therefore, the technique 
which could produce such a result (like the Borel summation of the "leading" asymptotics) can 
be quite misleading for calculating the low-energy spectrum as we know it from phenomenology. 

5 On the determination of the strange quark mass 

In case of the (us) current the mass of the strange quark should be taken into account. This was 
discussed recently in Refs. || 0, §]■ ^he conclusions we have previously drawn about moments 
with weight function (1 — s/M 2 ) k for the massless case remain the same. We comment only 
on a new feature of the massive correlators with the s quark mass m s considered as a small 
correction. The occurrence of the strange quark mass leads to an additional freedom in the 
choice of the scalar function for the momentum-type analysis because the correlator is no 
longer transverse (unlike in Eq. @). In the literature there are different suggestions about 
which of the correlator parts U q (q 2 ) and U g (q 2 ) in Eq. (|j) is better suited for the extraction 
of the strange quark mass. It is worth to have a closer look at the situation. While in present 
calculations radiative corrections are taken into account, the important term is still the leading 
order contribution. Because m s is small, this gives a correction which should be taken at the 
lowest order available in order to have a reasonable accuracy. The leading order spectral density 
for the current ifj^j^s related to the part H q (q 2 ) of the correlator in Eq. @ reads 

P 9 (s)={l-^ + ^jB{s-ml). (91) 

This spectral density is perturbatively insensitive to m 2 s because at large s the expression in 
Eq. ([91]) starts with a term proportional to m\. A term of order m 2 appears only after the 
integration of the spectrum over s with weight functions which are hard (nonvanishing) at the 
origin (at s ~ m 2 or s = for small m 2 ). Moments given by powers of s (as s l , I > 0) have no 
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contribution of order m 2 s . Indeed, for the zeroth order direct moment one finds 
where a term of order m 2 s occurs. On the other hand, for the first moment 

* - 2 € fe) P -W - 1 + » (I) 2 - 4 (I) 3 + 6 (S) 2 '" (I) (93) 

there is only a logarithmically enhanced contribution. For higher moments there are only 
higher order corrections, e.g. for I = 2 one has 

Therefore, perturbative s'-moments of U q (q 2 ) have no sensitivity to m 2 s . This is also expected 
for this correlator part because the mass corrections start from power corrections and lead to 
no cut which is a perturbative singularity in QCD and reliable for large s. Note in passing that 
this is also the reason why the perturbation theory expansion for the coefficient function of m 2 s 
in the correlator part U q (q 2 ) is known by one term less than that for the part U g (q 2 ): for the 
correlator part U g (q 2 ) only divergent parts of the four- loop diagrams are necessary while for 
the part U q (q 2 ) one has to calculate the finite contributions of the four- loop diagrams which is 
beyond the scope of present computational algorithms. 

The afore mentioned structure of the spectral density also means that there is no way to 
improve the theoretical accuracy of the contribution of the mass by making the corresponding 
coefficient function more perturbative: as soon as one tries to make the coefficient function 
multiplying m 2 more perturbative by introducing an additional suppression by a power of s (or 
by increasing I for the modified moments), the contribution of the mass becomes parametrically 
smaller, so instead of m 2 one has m\ in the leading order. When the anomalous dimension of 
the mass is taken into account, the power s then leads to an additional suppression by a power 
of a s . 

On the other hand, the (M 2 — s) k -moments all have the same contribution of order m 2 s given 
simply by the zero order direct moment Mq of Eq. (|92|). 

il; - w) k - 1 - 3( * + 4 + (S) • <"> 
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This amplitude is used in Ref. || where a better stability in the MS scheme was detected. 
In contrast to p q {s), for the spectrum p g (s) of the amplitude Ii g {q 2 ) in Eq. (0) one has 

/ \ / 3m 2 \ „ . 0n , . 

P a (s) =[- s + ^ ± -^) °( s ~ m ") ( 96 ) 

which is sensitive to m 2 for large s. Any perturbative moment of this amplitude has m 2 

contributions and can be used in the perturbation theory analysis of m 2 corrections. Namely, 

n ^ f M? f s V ^ ds n,2i* + l 3m 2 fm 2 a \\ 

M - {l - l) L {m?) p ^r- M n—2-^ +o {^))- (97) 

Therefore, the choice of the correlator part U g (q 2 ) is definitely more reliable perturbatively 
for analyzing mass corrections, a fact that is also discussed in Refs. ||, [?[]. However, the 
experimental situation is again less favourable to the choice of U g (q 2 ) in comparison to U q (q 2 ). 
The correlator part U q (q 2 ) contains a contribution of the K meson which is well pronounced 
and can be easily seen while the amplitude U g (q 2 ) contains only spin-one states and is poorly 
measured. In order to benefit from a better experimental accuracy, the use of the amplitude 
Il 9 (g 2 ) is favoured, but this choice leads to sum rules which are rather similar to the pseudoscalar 
sum rules (see e.g. Ref. ||54{| ) and does not provide new aspects for the strange quark mass 
determination specific to r decays. 

6 Integration of running effects to all orders 

In this section we discuss the technique of resumming the effects of renormalization group 
caused running of the strong coupling and other parameters (simply called "running" in the 
following) by taking the integration path along the contour of a circle in the complex q 2 - 
plane. As was discussed earlier in Sec. 3, the integration along the contour of a circle in 
the complex g 2 -plane is completely equivalent to the integration along the cut if finite-order 
perturbation theory expressions for the polarization functions are used. This is the direct 
mathematical consequence of the analytic properties of the functions ln p (— M 2 /q 2 ) for positive 
integer values p. However, if the renormalization group improved polarization function is used 
for the integration, some new features appear. This was also discussed in Sec. 4. The reason 
for this is that the analytic properties of resummed polarization functions are different from 
those in finite-order perturbation theory. 
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The integration along the contour including a full renormalization group resummation is 
now the most popular technique of accounting for the running of perturbative quantities: it 
efficiently resums an infinite number of terms generated by the evolution of the coupling con- 
stant ]55], [SfJ. Clearly it also includes some nonperturbative features when an infinite number 
of terms is resummed which means that this also constitutes a special recipe for resummation of 



perturbation series |57[. The formulae for the integration along the circle are known in a closed 
form as one-fold integrals but the final integration itself is usually done numerically. An analy- 
sis of the analytical form of these integral expressions is useful for clarifying the real content of 
such a resummation technique. In doing this we discuss how our general considerations given 
in the previous sections work within this resummation technique. It should definitely give the 
results confirming the physical considerations (effective scale structure) qualitatively, otherwise 
the technique of resummation would be irrelevant for physical applications. As we shall see, the 
resummation technique fully agrees with the conclusions about perturbative calculability and 
asymptotic structure of the moments found in the analysis on the cut in the previous sections. 

6.1 Quantities on the circle 

As in the case of the analysis on the physical cut, our main concern is to account for the 
running. Therefore, the introduction of an effective coupling is helpful. Here we introduce an 
effective coupling in the Euclidean domain which requires some new notation. 

We start with Adler's function in Eq. fljjD, introducing the effective strong coupling o:e(Q 2 ), 

d{Q 2 ) = -Q 2 ^n(Q 2 ) = 1 + (98) 

dQ Z 7T 



The effective coupling o:e(Q ) obeys the renormalization group equation (see e.g. |58| ) 

Q2 w (^P) = 0{aE(Q2)) - (99) 

For the considerations we are aiming at it is enough to use only leading order in the /3-function 
which contains already the bulk of the whole effect. Effects due to higher order corrections of 
the /3-function are really small and do not change the basic picture. They only slightly affect 
the conclusions numerically [[^J. Therefore, we consider the renormalization group equation 
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The resummed correlation function reads 



II(Q 2 ) = ]n(^)+i-ln( aE ^ Q2 n + subtractions (101) 
\Q J ft I M 



where 



Oi 

aE{Q) = i + w a T MHQVM ? ) (102) 

with a T = afi(M^). Taking this into account and parameterizing the contour by Q 2 = M 2 e lip 
one obtains 

Il(M 2 e ilp ) = -i<p - ^- ln(l + ip a T (p/n) + subtractions (103) 

Po 

where appropriate subtractions are added. The first term is the parton contribution which is 
independent of a T . This term is readily be taken into account. We therefore concentrate on 
the contribution coming from the second term which is a T dependent. 

The analysis of this second part can now be done for the moments Mm, given by 

M kl = l + mu = ^^^V./l 1 ^ / n(M*z)(l + z) k z l dz = 

2m kW. J\z\=i 

(k + I + 1)\ r 



r n(M r 2 e i¥, )(l + e i,fi ) k e iil+1)v, d<p. (104) 

J — 7T 



2tt k\l\ 

Taking the polarization function as given in Eq. (|103| ), the first part leads to the parton result 
Mki = 1 while the second, a T -dependent part gives the moment m^i where 

m kl = \> 1 h / (1 + e*)^ ln(l + ip a T cp/7r)e^dcp. (105) 



2tt/? k\U 

Note that the estimate of the saturation region for the moments given by integrals like the one 
in Eq. (|105|) is a bit trickier than in the case of the integration on the cut. On the cut the 
measure is positive and is normalized to 1 for any type of moments with n > 0, 

(n + 1) / x n dx = 1 and (n + 1) / (1 - x) n dx = 1 (106) 
Jo Jo 

which makes the estimate of integrals 

(n + 1) f 1 x n f{x)dx (107) 
Jo 

for continuous (and rather smooth) functions f(x) straightforward. In the present case, how- 
ever, the measure of the integration region is oscillating and vanishes for constant functions 
(k > 0, I > 0), 

(1 + e^fj^e^dy = 0. (108) 

J — 7T 
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Therefore, some additional care in the consideration of essential regions of integration is nec- 
essary for the formulation of moments on the contour. Note that we have already encountered 
such a problem while analyzing the moments of the narrow resonance contributions. Therefore, 
the intuition gained from the situation with continuous functions generally does not work for 
distributions. 



6.2 Techniques for the calculation of contour moments 



Because explicit functions are given in Eqs. (|101|) , ( p.03| ), (|104|) , and ( |105| ), the discussion of the 



representation and properties of the moments on the contour is a pure mathematical problem. 
Since that moments on the circle (as in Eq. ( |105[ )) are just Fourier coefficients of the correlation 
functions U(Q 2 ), the technique for working with these is well-known. One finds similarities 
to the evaluation of sunset diagrams for the case when the problem is reduced to a Fourier 
transform. 

For /?o«t < 1 the moments in Eq. ( |105|) can be expanded in a convergent series in a T . The 
finite radius of convergence within this technique of resummmation is a general feature which 



persists in higher orders of the /3-function and was studied in some detail |59[| . The convergence 
radius decreases when higher orders of the /3-function are included. If we compare the value of 
the effective coupling in r decays (not to be mistaken with the coupling within the MS-scheme) 

0.14 (109) 



7T 

with 

-L = i- = 0.1415 ... , (110) 
irp 97r 

the relation ^a T < 1 is still valid, though almost marginally. However, the exact expression 
given in Eq. ( |105| ) without expansion of the logarithm provides an analytic continuation beyond 
the convergence radius even when a T lies outside the convergence radius. 

We first consider the basic moment moo- The discussion of this basic moment leads to the 
technically simplest (shortest) expressions while containing already all features of the general 
case. The generalization to other moments (and types of moments) is then straightforward. 



The most direct way of calculation is to expand the logarithm in the integrand in Eq. ( |105|) 
in a series in a T which gives nothing new in comparison with the finite-order perturbation 
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theory case. However, one can generate all terms of the series because the function is known. 
The series obtained by this procedure represents an analytic function of the coupling constant 
at the origin. Therefore, for small a T this series can be used for a numerical evaluation of the 
moments. This procedure may not be the best one concerning economy of evaluation since many 
terms may be needed to determine values of a T close to the convergence boundary. A finite 
piece of the convergent series gives different accuracy and can also completely misrepresent the 
function if a T is beyond the convergence radius (which is not so dramatic for positive values of 
a T , though). 

One can proceed with the analysis of the moments in a different way by constructing just 
an efficient computational scheme. Integrating n times by parts one obtains (we remind the 
reader that we use the zeroth order moment for simplicity) 

m m = —\<p+}_^{]-l)\\ smj0+ / -—Mill 

TTPO { ~[ \ Tir J 2 \ 7T ) J-n {l + i(3 a T (p/iT) n ) 

with the polar coordinate functions r and 4> (note the difference to (p\) defined by 



1 ± i(3 Q a T = re ±l(t} , r = ^/l + ffia 2 T , <f) = arctan(/3 a r ). (112) 

The n-fold integration by parts removes a polynomial of order n from the expansion of the 
logarithm in Eq. ( |103| ). 

One can see that the result is an asymptotic expansion as the residual term, i.e. the last term 
in Eq. (|111|) is of the formal order a™. However, the obtained result is not a series expansion 



in the original coupling determined in the Euclidean domain but a more complicated system 
of functions related to it. The system of functions is ordered and the asymptotic expansion 
is valid in the sense of Poincare. The system of functions is obtained by using the expression 
for the running coupling in Eq. (|102|) in the Euclidean domain and continuing it into the 
complex plane and onto the cut. When the analytic structure of the initial function is known, 
asymptotic expansions which converge fast for the first few terms (as a representation in the 



form of Eq. ( |1 1 1| ) ) are more useful for practical calculations than formal convergent series that 
require many terms for getting a reasonable accuracy. 

The expansion in Eq. ( |1 1 1| ) can give a better accuracy (for some n and a T ) than a direct 



expansion in a s . Indeed, this expansion includes a partial resummation of the 7r 2 terms which 
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is a consequence of the analytic continuation p9j. Therefore, the expansion can be understood 



as being done in terms of quantities defined on the cut. This is evident because in some sense 
the procedure consists in calculating derivatives of the spectral density on the cut. Because the 
region near the real axis is important, the continuation causes a change of the effective expansion 



parameter a T — > a T /y 1 + Po a r- The first term in the expansion shown in Eq. (|111|) is just the 
value for the spectral density expressed through the coupling in the Euclidean domain. It can 
be considered as a change of the scheme at large values of s. 

With the concise expression for the moments at hand one can change the form of the residual 
term. The relation 

'{3 a T \ n r* e^dcp 



n-l)\ 



7T 



-7T (1 + if3 a T ip/ir) n 



7r J \J-oo Jw J (1 + i/3 a T (f/n) n 
valid for any n leads to a representation of the zeroth order moment in the form 

m 00 = ^-l^^ + ^ + EO'-l)^— YsinO^) 
vr/3 1 £i \ nr J 

(n-1)! (Po<* T Y{ r * + r°°\ tPdip \ 



2 y 7r y \7-oo A / (1 + i(3oa T ip/iry 
Here an explicit nonperturbative term e~ n ^° aT has appeared. Eq. (|114 ) and Eq. ( p. 11] ) are 



formally different but actually identical. Therefore, the choice for the expansion (or repre- 
sentation) for the moment is a question of calculating the residual term. This quantitative 
consideration is, of course, only possible if one has a concise expression for the function from 
which the series is generated as given by Eq. (|105|) in our case. Any conclusions about the 
precision or the analytic structure of the sum of the series based on the terms of the series only 
without a specification of the residual term are rather useless as one can see from Eq. (|113|) . 

We stress that the moments in Eq. ( |105[ ) are analytic functions of a r for small values of the 
coupling a T . This means that the non- analytic piece in Eq. (|113| ) cancels the corresponding piece 
in the residual term. If the residual term is dropped, the analytic structure drastically changes 
depending on which representation, either Eq. ( |1 1 1| ) or ( |114| ), is used. This demonstrates the 
danger of reaching conclusions about power corrections emerging from the extrapolation of the 
running to the IR region. Mathematically, or from a theoretical point of view, a conclusion of 
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Figure 7: the integration contour in the complex yj-plane 



this type is a hypothesis, especially having in mind that such power corrections are not necessary 
in any way for the consistency of the theory (in contrast to the OPE corrections which are well- 
defined in perturbation theory and are motivated by, say, explicitly distinguishing different 
channels) . Physically, the use of such corrections is very difficult to appreciate if higher order 
terms of the expansion in a s are taken into account, i.e. an explicit nonperturbative term in the 
expression, not being special in any general way, does not numerically differ from the higher 
order terms as well. 

Yet another representation for the moments can be readily obtained. It is possible to recover 
the form of the moments as integrals over a spectral density (see below). But because moments 
are calculated as explicit functions, we can obtain a result by using a simple change of variables 
in a pure mathematical sense. In order to obtain this new representation we go to the complex 
plane in (p (see Fig. 0). As a technical trick we first integrate once by parts to obtain 

a T r (1 + e^)dip 
m 00 = — — — — 115 

2tt z J --k 1 + ipoa T tp/ii 



This representation is somewhat different from the representation in Eq. (|1 1 1|) for n = 1. They 
differ by an integral which can be explicitly computed, 

^ ■ In ( t— '.'f ° ) = — 5- arctan(/? «)- (H6) 



2it 2 J -it 1 + if3 a T (p/n 2mf3 \1 — if3 a a T ) %/3 

Now we consider the integration over a rectangular contour in the complex <^-plane. The 
part of the contour on the real axis from — 7r to 7r leads to the moments. The integral over the 
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contour is given by the residue at the pole <p = i7i//3 a T . We thus have 

a T r (l + e itp )d(p 

m 00 



+ 

ip=in//3oa T 



2tt 2 J~n 1 + i(3 a T (p/ir 
ia T ( 1 + e ilp 
7r \i(3 a T /ir 

a T f°° (1 + e™-t)id£ 



roc [i + e-^m , ,- n 

2tt 2 Jo 1 + i/3 a T - /3 a T Z/ir [iP ^ 

Z7T Z JO 1 — «Po«r — PoCtr^/ll 



1 + e^- R )d(p 



+ 2^tLL l + i( 3 a T ^- ( 3 a T R/n ^ = * + ^ (U7) 
where we have indicated the substitutions for the different parts of the path. The last part 
vanishes because the integrand vanishes for R — > oo while the integration range is finite. We 
thus obtain 

id + e -/*«o 4f n (1 ,~ e ~*?5 ^ (us) 



m °° /V~ ' " ' /Wo 7r 2 + (£-7r/A,a r ) 2 ' 

Next we replace —7[//3 a T = ln(A 2 /M 2 ) and substitute — £ = ln(s/M 2 ), so that = s/M 2 

and — <i£ = ds/s. One obtains 

A further transformation gives 

1 /A 2 \ 1 f M? ( ln( S /A 2 ) \ ds , 1on . 

moo = _ — - -\ / arccos v ' ' — -. (120) 

fo\M*) nP Jo {^ 2 + \n 2 (s/A 2 )) W 1 ' 

One easily recognizes this representation as an integration over the singularities of II(g 2 ) in 
Eq. (|101|) . In addition to a cut along the positive semi-axis there appears also a part of the 
singularity on the negative real s-axis. This part is a pure mathematical feature of the concrete 
approximation chosen for II (q 2 ) and is not related to the physical content of the problem. 
Indeed, if moments are written as explicit functions we can calculate them in the way we find 
most convenient for a concrete application. The result reads 

m oo = / , f , (121) 



with 



Ml 



a(s) = -^(A 2 + s)6(-s) + JU( S ) arccos f ln(g/ f 2) | (122) 

A) tt/5o V\/ 7r2 + ln ( S M 2 ) 
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This formal result can be reformulated as integration over the spectrum a(s) using Cauchy's 
theorem. Indeed, from Eq. ( |101|) one readily deduces the singularity of the resummed polariza- 
tion function with the discontinuity 

Disc nO) = ^ \e(A 2 + s)9(-s) + -6(s) arccos f M^/A 2 ) _ \ 1 _ 



which coincides with a(s) in Eq. ( |122| ) and, thereby, simply represents the spectrum. 

This spectral function a(s) is shown in Fig. |8|. The part of the spectrum on the positive 
real axis is an analytic continuation of the renormalization group improved correlator function 
Tl(Q 2 ) to the cut |3^, [55|, 60, 61, |62 |. It can be conveniently written in the form 



, ln(s/A 2 ) \ . / vr \ / vr \ 

arccos . = arcsm . = arctan — — — — — - (124) 

/7r 2 + ln 2 (s/A 2 )/ Vv /7r2 + ln2 ( s / A2 )/ V ln ( s / A )/ 



with the appropiate branches of the relevant functions. Introducing 



one can write 



a(s) = — 7~7~i o\ (125) 
/3 ln(s/A 2 ) 



arctan ( ) = arctan(/5 o;(s)) (126) 



V ln(s/A 2 ;/ 

Taking the value of the spectral density at M 2 as initial value 

1 in \ a r , s^t 2\ 



a(M 2 ) = — arctan(/? a T ) = — + 0(a 2 ). (127) 
7T/? vr 

we can introduce the continuum part of the spectral density, 

a c( s ) — —pr arctan(/3 o;(s)). (128) 

7T/3 

The differential equation determining the continuum part a c (s) through its initial value a(M 2 ) 
can be constructed by differentiating Eq. (|128|) with respect to s, 



as y 7T y 1 + pQa(s) 2 

with a(s) taken from Eq. ( |1 25|) . By inverting Eq. ( |128| ) we have 

Poa(s) = tan(7r/3ocr c (s)) for s > (130) 
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Figure 8: The spectral function a(s) as a function of s. Shown is the exact result for the 
discontinuity (full line) as well as results for non-closed circle integrals approaching the axis 
along a full circle from both sides, see Fig. 0. The numerical values for the distances orthogonal 
to the real axis towards the positive and negative imaginary semi-plane in s are given in units 
of Ml 



and therefore obtain 

s-^-cr c (s) = }— sm 2 (nP a c (s)) for s > 0. (131) 

as it Po 

This equation can indeed be considered as an evolution equation for the spectral density a c (s) 

determining er c (s) through its initial value cr(M^). Therefore, one can simply introduce an 
effective charge 

oim(s) = 7ra c (s) or a M (s) = — ^— = a c (s) (132) 

7T 



(see Eq. flUf)). The evolution equation for a M (s) is then given by 

sm 2 (wp a M (s)). (133) 



da M (s) 1 . 2 



ds K 2 f3 

Therefore, we now define the coupling as the value of the spectral density on the cut far from 
the IR region. This is a perturbation theory definition. The evolution of this coupling, however, 
is calculated by taking into account the analytic continuation. Then it has an IR fixed point 
with the coupling value 

a M (0) = I/fa (134) 
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(we assume that all transformations are made for small au so that this is the first IR point). 

Starting from finite-order perturbation theory one therefore can find the all-order result 
because of analytic continuation and running in the IR domain - it happens to be finite (no 
poles for the spectrum at small values of s). Therefore, the function a c (s) can be used to 
generate moments by integration on the cut. Note that in this approach the part of the 
spectrum at negative s is lost and cannot be recovered from finite-order perturbation theory. 

Note, however, that the value at zero mentioned in Eq. ( |134j ) is not a universal quantity. If 
Adler's function starts with a proper power of the coupling constant as it is the case for gluonic 
observables, for instance, this picture will change. For instance, we take 



DiQ i ) = (e^m\\ (135) 



Then the corresponding spectral density in leading order /3-function approximation reads 

1 1 a 2 (s) 

P{S) = M hi 2 (s/A 2 )+7r2 = 7T 2 (l + /3 2 a 2 (s))- (136) 

As in previous cases, we can define an effective coupling in the Euclidean domain, 

a M (s) = . a{s) (137) 
n^l+P 2 a 2 (s) 

For this effective coupling we indeed have (cf Eq. ( |134|) ) 

a M (0) = 0. (138) 
The /3-function for the effective coupling obtained from Eq. ( |137| ), 

Pm^m) = -f3 a 2 M ^l - (n/3 a M ) 2 = ~Poa 2 M (l - ^(7r(3 Q a M ) 2 + O(o^)^ (139) 
differs from the /3-function of the effective coupling clm (cf. Eq. (|133|) ) 

l3 M {a M ) = sin {irf3 a M ) = -(3 a 2 M (l - \(-K^a M ) 2 + 0(a 4 M )) (140) 

7T z po V 3 / 

at next-to-leading order. 

The consideration can be generalized to higher powers of the coupling in the Euclidean 
domain. For 

Dpm = (141) 
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Figure 9: The spectral function p p (s) as a function of s, given by Eq. (|142 ) . All spectral 
functions for p > are going to zero for s —>■ 0, for p > 2 they have nontrivial zeros and a 
fluctuating behaviour near the origin. 



we obtain 



where 



1 sin p0 s 

M s = — n 142 



r(s) = /3 Vln 2 (s/A 2 ) + tt 2 , tan0(s) = ^ ^ = /%a(s). (143) 



For p = we retain Eq. (|128 



Po( s ) = —r^i 3 ) = —T arctan(/5 a(s)). (144) 
7T/?o 7t/? 

But for large values of p the function p p (s) starts to fluctuate when starting at p(M 2 ) and 
approaching s = where it finally reaches zero. Because if this behaviour, shown in Fig. |9] in 
some detail for the first four values of p, p p (s) fails to be interpreted as a (positive) spectral 
density in a region which enlarges with increasing values of p. 

We conclude that the two recipes of resummation (on the contour and on the positive 
semi-axis for s) are different: they differ by the integral over the negative real semi-axis for 
s or, more generally, by the integral passing through the infrared region. This point is worth 
discussing in more detail. In the contour formulation it is not essential what particular point- 
by-point behaviour in the IR region exists. For the analytically continued correlator this is not 
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Figure 10: Contours in the complex plane with a fixation in the region A, taking into account 
possible occurrences of singularities. The figure part on the left hand side shows the standard 
circle path which circumvents the singular region B while the occurence of other singularity 
regions as discussed in the text (regions B' and B") may lead to different possibilities for 
choosing a path (C a , Cb, C c ). The path Cd crosses the singularity region and, therefore, cannot 
be used from the perturbation theory point of view. 



important unless the contour crosses a nonanalytic region. Whatever singularities exist in the 
IR region (Regions B, B', or B" in Fig. [TO), the contour includes them. This is the definition 
of resummation on the contour which is explicitly perturbative. 

For the formulation on the cut the extrapolation of the running of the coupling constant to 
the IR region is crucial because it provides an interpretation of the spectral density in the IR 
region. One can formally write 

a(s)ds (145) 



and try to interpret this integral. The naive substitution of the expression in Eq. (|125|) leads to 



a singularity which is known as Landau pole. Still, further formal manipulations readily give 
(t = /5 a(M r 2 )ln(M T 2 / S )) 

Ml r M * q(M?)ds _ M 2 e^'^^dt 

a[S) Jo l + /3 a(M 2 )ln( S /M 2 ) " ATio ~t " (M6) 

Both equations are not well-defined from the beginning, the second form being a Borel rep- 
resentation. The same problem can be reformulated as a divergence of the asymptotic series. 
Indeed, by expanding the expression for the running coupling under the integration sign in a 
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perturbation theory series one has 

The summation of the series in Eq. (|147Q is equivalent to the interpretation of the integral, 
independent on what form is chosen for it, so either Eq. ( |145| ) or Eq. (HE). Therefore, an 



integrable behaviour of the coupling constant at small s (whatever the definition of this quantity 
is) can solve the problem of divergence of the integral or the summation of the asymptotic series. 
But this solution is strongly model dependent because the extrapolation of the evolution into the 
IR region is essentially arbitrary and higher order corrections are important. The explicit form 
of the extrapolation in Eq. ( |122| ) gives an extrapolation motivated by analytic continuation. It 
can also be considered as a special change of the renormalization scheme [p6 |. Indeed, for the 



coupling au with the evolution given in Eq. (|133|) one obtains an IR fixed point. In any case 
the result depends on the IR region which is strongly nonperturbative. 

As it has become popular as being used as an extrapolation, for the analytic continuation 
the moments are defined simply by 

anal °c(s)ds 
m 00 - /_ T7Z , I 148 ) 



M r 2 

excluding totally the negative part of the spectrum that emerges in the exact treatment of 
quantities originally defined in the Euclidean domain (see Eq. ( |121|) ). Using Eq. ( |120|) one can 
relate these moments to the moments on the contour by 

m ^ al = e -*/A>«r + moQ , ( 149 ) 

The contour definition for the moments is a strictly perturbative one and in this particular 
case simply gives the analytic function of a T at the origin a T = 0, i.e. moo(« T ) is analytic at small 
a T . Then the decomposition in Eq. (|149 ) is at least meaningful, the right hand side does contain 



terms of different analytic structure. We can, of course, reproduce the contour moments by 
summing on the cut (which gives m^ 1 ) and adding a pole (giving the "nonperturbative" term 
in Eq. ( |149| )) which is an allowed extrapolation as well (because the IR region behaviour cannot 
be found in perturbation theory). But it is difficult to imagine that one could suggest such a 
treatment of the situation just as it is, if not derived naturally from the contour definition. Note 
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that this is completely analogous to the situation with Coulomb resummation: by summing on 
the cut 4m 2 < s < oo one can never restore the pole contributions generated by the bound 
states. Only an explicit analytic continuation can reveal the whole analytic structure of the 
correlator. The result of Eq. ( |149|) can readily be generalized to any power of basic 
quantity which is important for correlators of hadronic currents with nonvanishing anomalous 
dimensions. 

Note that the problem of a particular choice for the extrapolation of the evolution into the 
IR region cannot be uniquely solved from the point of view of the perturbation theory because 
the IR region is nonperturbative and there is no tool to choose a "correct" extrapolation in 
cases when there is no explicit analytic structure known. Therefore, the extrapolation is rather 
arbitrary for basically all cases except the case of two-point correlators. Still, a particular 
choice for the extrapolation can be important numerically for the theoretical description of 
the available high precision data. In r decays the special role of contour moments is related to 
perturbation theory - they are really computable in perturbation theory without any additional 
hypotheses. For the s'-moments the problem is less important numerically while for the (M 2 — 
s) fc -moments it is crucial (see Table [|). This is obvious because the (M 2 — s) fc -moments are 
saturated at small s and the particular way of treating this region is important. In order to 
discuss this issue, we now leave the special case k = I = and turn to the description of 
modified moments, discussing the differences. 

6.3 The moments m$i 

In general, for any given moment m ; one finds 



The representation obtained by using integration by parts shows an improvement in the con- 
vergence for moments at large I equivalent to the replacement a T — > a T j il + 1). This coincides 
with conclusions drawn from finite-order perturbation theory. For arbitrary I one can also find 




(150) 
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a representation containing an explicit nonperturbative term, 

7T/3 { ~i \(l + l)irrj 

l+1 (n-l)\ [ /3 a \ n ( f~* , /-°°\ e'C+^tip 



+(; + l)(-ir^ & j ^ + y j__r_|. (151) 

These formulae give a direct generalization to large I moments and can be treated perturbatively. 



6.4 The moments rriko 

The moments for I = corresponding to the a T dependent part of the correlator in the leading 
order are given by 

m k0 = f (1 + e^) fc ln(l + i(3 a T <f/7r)e^d<f. (152) 

These moments can also be expanded in a convergent series in a T for small values of oe T , 
recovering the results from finite-order perturbation theory. 

The usual arguments given for the preferred use of these moments from the theoretical point 
of view is that they suppress contributions close to the real axis where OPE is not applicable 
(region A in Fig. [TDD). They do suppress these contributions because of e lip ~ — 1 in this region. 
One can also obtain this result by looking at the integral measure for the moments m^o, 

(1 + e^e^dip = (2 cos(y>/2)) V fe ^V<^ (153) 

where the suppression is given by the factor (1 + e lip ) k . The suppression of oscillations starts at 
larger values of the integration parameter (p than for the moments rrioi but there is a suppression 
by the power of cos(<£>/2), cos(<£>/2) = at ip = ±7r. However, one can see that these moments 
are almost insensitive to the value of a T for large k (or, more precisely, to the perturbative 
structure). It is clear that these functions determine integrals which really have no relation 
to perturbative running because such a running is completely suppressed. This makes an 
asymptotic expansion in a T almost useless (we could not find any efficient asymptotic expansion 
in a T for this moments at large values of k). However, it is straightforward to write down a 
nonperturbative representation for these moments at large values of k. A representation in 
terms of integral over the spectrum is particularly convenient. One has 



m k0 



(1 + e^) k ln(l + H3 oaT <p/*)e t *d V > = / . „ 1 - ^ (154) 



2vr/5 J-* ' K ^ ^' ' r J-A* V J M r 2 
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k = 1 


k = 2 


k = 3 
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k = 6 


k = 7 



Table 5: Contribution of the integral taken over the interval from —A 2 to +A 2 for the moments 
rriki in Eq. ( |159| ) relative to the integral taken over to the whole integration range. 



The part of the spectrum on the negative real axis is given by a constant. For large values of 
k the integral is saturated at the left boundary of the integration region and reads 

/o 



1 



1 







Ml 



1 +(*+i) r(i 



s \ a c (s)ds 
M 2 / 



Ml 



(155) 



The contribution of the integral is suppressed because the spectrum a c (s) is a smooth and finite 
function leading to a finite result for the integral. 

Another way to look at the problem is given by integrating by parts once. One has 



1 



dip 



27r 2 1 + i/3 a T (p/iY 
In this case the boundary terms vanish. One obtains 

1 r M r 



(156) 



1 / A 2 \ fc+1 
mk0 = Jo{ l + K) 



1 



\ k+l 



ds 



f3 Jo \~ M^j 7T 2 + ln 2 (s/A 2 ) s ' 
Here it is easy to estimate the contribution of the integral as well. One has 



(157) 



k+l 



1 - 



ds 



< 



Mi 



ds 



< 1. 



(158) 



M 2 J 7r 2 + ln 2 (s/A 2 ) s Jo vr 2 + ln 2 (s/A 2 ) s 
Note that the whole function rriko is still represented by a convergent series for small values of 



In general we determine the contribution of the nonperturbative region determined by the 
scale A 2 to the modified moments 



[ M * (1 _ jj\ k ( s V a ^ ds 

1-v \ Ml ) \M1) Ml ' 



(159) 
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Table 6: Contribution of the integral taken over the interval from —A 2 to for the moments 
rriki in Eq. ( |159| ) relative to the integral taken over to the whole integration range. 



in order to consider the enhancement of this region in relation to the moment as a whole. In 
Table |5| we give values for the contribution of the interval range [—A 2 , A 2 ] to the contribution 
for the whole range [—A 2 , M 2 ] for different values of k and I, taking A 2 = 0.2 GeV 2 . A strong 
enhancement can be observed especially for the case / = and large values for k. One observes 
that the first term is larger than the second one already for k > 3. Though weakened, the same 
tendency can be seen in Table |6] where we considered the contribution of the interval range 
[—A 2 , 0] relative to the whole range. 

This observation also suggests the question of sensitivity of the moments to the parameters 
describing the strength of strong interactions, namely the values of a T M and A. As expected, the 
s'-moments are sensitive to the spectrum at large energies, i.e. cx(M 2 ) which is almost equal to 
a T M . The large k moments are directly sensitive to A as one can see from Eq. (|157|) . However, 
perturbation theory is insufficient for these moments, and thus the theoretical formulae for large 
k moments like Eq. (|157|) cannot be confronted with experiment and should be supplemented 



by a quantitative hadronization mechanism. 



6.5 Higher orders of running 



The formulae obtained for the polarization function can readily be generalized to higher orders 
of the /3-function. For Adler's function one has 

D{Q 2 ) = l + a E (Q 2 ). (160) 
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The polarization function is then determined by the equation 



d 
dQ 2 



-Q 2 -^i{Q 2 ) = l + a E {Q 2 ). (161) 



A simple way to calculate Tl(Q ) is to use the ansatz 

n(Q 2 )=ln(j^) + \nf(a E (Q 2 )). (162) 
One then obtains a differential equation for f(a E ) which reads 

- Q2 >^ = 1 - n^mfw ! ^ = 1 - 7^w%^- < 163 > 

Using the renormalization group equation 

Q^-^f 1 = P (MQ 2 )) (164) 

one finds 

/' (a E (Q 2 )) P (a E {Q 2 )) = -f (a E {Q 2 )) a E (Q 2 ). (165) 
This differential equation can be solved in quadrature, 

"^-/tS- (166) 



The solution for H(Q 2 ) is therefore given by the trajectory 

a E d(L E 



(167) 



/3(a E ) ' 

Eq. (|167|) can be used for any f3- function (with the restriction that no singularity should occur). 
For /3(a) = —floa 2 we obtain the former result, 

(3{a) = -(3 a 2 \n f = (a E (Q 2 )) . (168) 

For P(a) = —/3oa 2 — (3\a z we explicitly obtain 

-«n-(SH h (s5s&j)- 

The generalization to any polynomial function that may occur in perturbation theory is straight- 
forward. 



64 



As a last note we comment on the Landau singularity which occurs at the leading order 
of the running. For two-point correlators the problem is completely solved by the contour 
integration to any order of the /3-function. We stress that the pole is inside the circle (see 
the discussion in Ref. |57]]). One cannot deform the contour to be close to the origin because 
this region is completely nonperturbative and should therefore be avoided. This situation is 
completely analogous to the situation with Coulombic poles [^3[ or with a spectrum below 
threshold for a gluonic operator induced by a heavy quark loop |M| . Any type of singularity is 
avoided by moving the contour far from the origin and keeping the singularity inside, thereby 
including also the integral of the spectrum. 

7 Conclusion 

We have shown that the direct s'-moments are perturbative and can be reliably calculated 
within operator product expansion (OPE). The mixed moments with the weight functions 
(1 — s/M 2 ) k are strongly nonperturbative. This is natural and is expected in QCD with asymp- 
totic freedom. We have analyzed finite-order perturbation theory, power corrections and re- 
summation. The results agree with each other and the picture is completely consistent. An 
asymptotic analysis (large values for / or k) definitely disfavours the large k moments despite 
arguments that they suppress the non-OPE region along the contour - these arguments are 
not valid. In practical applications to r decays and the determination of m s the second mixed 
moment with weight function (1 — s/M 2 ) 2 is sometimes used. Along with the initial phase space 
weight (1 — s/M 2 ) 2 (see Eq. (|])) the total results gives a moment with the weight (1 — s/M 2 ) 4 
which is marginal in the sense that it enhances nonperturbative contributions. Our consider- 
ation show that this moment is really on the border of what is still theoretically computable. 
The theoretical error becomes larger than one would expect for such an analysis. Note that 
an estimate of the systematic theoretical error, i.e. the error due to the functional dependence 
like perturbation series truncation, is complicated as compared to a simple statistical error of 
theoretical formulae due to uncertainties in the parameters used. The theoretical error for the 
moment with weight (1 — s/M 2 ) 4 appearing in the phenomeno logical analysis is big enough to 
be concerned. The point is whether we should worry about the size of this error in compar- 
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ison with the gain due to the improvement in the accuracy of the experimental moments or 
to trade this moment for a more reliable theoretical expression confronted with a less precise 
experimental moment. This still remains to be seen. 
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